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Several models of a strongly interacting Bose gas in an optical lattice are studied within the functional- 
integral approach. The one-dimensional Bose gas is briefly discussed. Then the Bose-Einstein condensate 
and the Mott insulator of a three-dimensional Bose gas are described in mean-field approximation, and the 
corresponding phase diagrams are evaluated. Other characteristic quantities, like the spectrum of quasiparti- 
cle excitations and the static structure factor, are obtained from Gaussian fluctuations around the mean-field 
solutions. We discuss the role of quantum and thermal fluctuations, and determine the behavior of physi- 
cal quantities in terms of density and temperature of the Bose gas. In particular, we study the dilute limit, 
where the mean-field equation becomes the Gross-Pitaevskii equation. This allows us to extend the Gross- 
Pitaevskii equation to the dense regime by introducing renormalized parameters in the latter. 
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1 Introduction 

The quantum statistics of non-interacting particles was established by S. N. Bose in 1924 QJ. Bose was 
able to deduce Planck's radiation law on the assumption that each quantum state can be occupied by an 
arbitrary number of indistinguishable photons. By applying this idea to the quantum statistics of an ideal 
gas of iVtot atoms enclosed in a volume V, A. Einstein predicted the occurrence of a phase transition 12: 
Below a critical temperature T c , a certain fraction of atoms would "condense" in the ground state of the 
system. This phenomenon is called Bose-Einstein condensation (BEC). 

In a homogeneous ideal Bose-gas (i.e., in the absence of an external potential), the critical temperature 
of the ideal Bose gas is given as 0] EJ |6] 0] 



where fee is Boltzmann's constant, h is the reduced Planck's constant, rttot = N/V is the particle density, 
m is the mass of the particles, and ((x) is Riemann's Zeta-Function. The condensate fraction is given as 



where hq is the condensate density. 

Historically, the first candidate for a possible realization of Bose-Einstein condensation was superfluid 
4 He, discovered by P. L. Kapitza in 1934 below T c = 2.2K. Although superfluid Helium is far away 
from the ideal Bose gas considered by Einstein because of strong interactions between the Helium atoms, 
the phenomena of superfluidity and BEC are related. Superfluidity was first explained by L. D. Landau 
in 1941 by an argument which is based on the idea that the viscosity of a fluid depends on the existence 
of quasiparticle excitations. Those excitations are created by friction between the fluid and a wall of the 
container. When the fluid has a velocity v relative to the wall, these excitations are relevant only if their 
creation at momentum k is energetically profitable, i. e. if the excitation energy is negative [4j: 

E k + hk ■ v < . 

Here Ey is the quasiparticle spectrum. In other words, the superfluid is destroyed by excitations if the 
velocity |v| exceeds a critical value v c with 



where the minimum is calculated over all the values of k. If the spectrum is linear for small momenta, a 
non-zero value of v c is found. It is important to notice that superfluidity and BEC are not identical. For 
instance, an ideal Bose gas can condense, but it is not superfluid due to Landau's principle, because the 
excitation spectrum is quadratic in k and therefore v c is zero. On the other hand, a weakly-interacting 
two-dimensional Bose gas satisfies Landau's criterion for superfluidity, but long-range order cannot appear 
due to the Mermin- Wagner theorem J8] [9j [TO], therefore there is no BEC. 

In an interacting Bose gas of uncharged atoms, the main contribution to the interparticle interaction 
comes from s-wave scattering between two particles. The characteristic length scale here is the scattering 
length a s . We assume a s to be positive, although it can also be negative in trapped Bose gases (without 
trapping potential a Bose gas with negative a s is instable H). For theoretical description, usually two-body 
interaction is assumed. Approximately, the two-body interaction potential can be written in the form of a 
(5-potential: 



2 




(1) 




(2) 



v c = mink — 



Mnt(r-r')«p<5(r-r') . 



(3) 
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Here, g is the strength of the repulsive interaction between two bosons. It is connected to the s-wave 
scattering length by the relation J4) 

g = ^JL. (4) 

TO 

This approximation is justified if the a s is small compared to the thermal de Broglie wavelength, the inter- 
particle spacing, and the characteristic length scale of the trapping potential [5]. It is possible to tune the 
scattering length over a large range of values (positive as well as negative) to reach the strongly interacting 
regime, where Bogoliubov theory is not applicable anymore ll4l [TT1 [T2ll . These magnetic Feshbach reso- 
nances became possible after the development of optical trapping as an alternative to magnetic trapping. 

After the introduction of an external potential V cx t, the full Hamiltonian of the Bose system in terms of 
bosonic field operators is 



H = d 



13, 



^ + (r) ( -^V 2 + V cxt (r)) 4>(r) + 9 - ^+(r)^+(r)^(r)V-(r) 



2to y / 2 



(5) 



The ground state of this interacting many-body system is not known, therefore the condensate density can- 
not be defined by the population density of the ground state like in the ideal Bose gas. An appropriate 
definition for a homogeneous system is the concept of "off-diagonal long-range order" which was devel- 
oped in the 1950's l4ll5l [T3ll . The condensate density is given by the long-range behavior of the one-particle 
correlation function 

n := lim (<$+ . (6) 

r — r' — >oo 

If the one-particle correlation function decays exponentially or algebraically, the condensate density is zero. 
An algebraic decay is found in a two-dimensional Bose gas at low temperature and in a one-dimensional 
Bose gas at zero temperature [ 14 1. 



1.1 Dilute Bose gas 

When the mean distance between atoms is large compared to their scattering length, which is the case when 
^tot^s "C 1, the system is said to be in the dilute regime. In this case, the effect of interaction is small. A 
consistent mean-field theory of a dilute Bose gas which is valid for low temperatures T<T C was given by 
N. N. Bogoliubov in 1947 ||3]|4]. The condensed phase is described by replacing the bosonic field-operators 
by the sum of a complex condensate order parameter $o an d fluctuations out of the condensate as 

$(r,i)=# (r,t)+^(r,i) , (7) 

where the field operators ip of the fluctuations fulfill bosonic commutation relations. This theory gives 
elementary excitations out of the condensate which have the energy spectrum 

where k is the wave vector. It is linear for small momenta ("phonon spectrum") and therefore satisfies 
Landau's criterion for superfluidity, in contrast to Einstein's non-interacting Bose gas with a quadratic 
energy spectrum. An important feature of an interacting Bose gases is the ground state depletion, which 
means that even at T = the condensate fraction is smaller than 1. This is also found in Bogoliubov 
theory. In a dilute Bose gas, the condensate depletion is small. 

The condensate order parameter $o 1S connected to the breaking of the global U(l) symmetry, which 
reflects the fact that the replacement 



(9) 
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where a is a global phase, does not change the physics of the system. The phase a can be chosen arbitrarily, 
but once it has been chosen, the symmetry is broken. This is the case in the BEC phase. This phase a is 
responsible for the fact that the quasiparticle spectrum in Eq. (0 vanishes for k = 0: The Goldstone- 
theorem states that the existence of a broken U(l) phase symmetry leads to a gapless excitation spectrum 

03. 

The order parameter is interpreted as a macroscopic wave function and can be split into its modulus and 
phase: 

$ (r,<) = |$o(r,t)|e ie ^) . (10) 
The local condensate density is related to the modulus squared of the order parameter 

n (r,t) = |$oM)| 2 , (11) 

and the gradient of its phase, V6*(r, t), is associated with the velocity field of the condensed atoms. Gross 
and Pitaevskii have independently derived an equation to describe the dynamics of the order parameter, 
which is known as the Gross-Pitaevskii (GP) equation |0][6][5]: 

-|^V 2 + F ext (r) + 5 |$ (r,i)| 2 ) $oM) = iftjUo(r,i) . (12) 

The third order term in &o, which is proportional to the interaction constant g, can be interpreted as the 
coupling of the order parameter to the local particle density as given in Eq. (fTTT >. For stationary solutions 
of the GP equation we use the ansatz $o( r , t) = $o( r ) exp(—i/j,t/h), where fi is the chemical potential. 
The GP equation then reduces to the stationary form 

l^-V 2 + K xt (r) -n + 5 |$o(r)| 2 ) $o(r) = . (13) 
1 .2 Trapped Bose gas 

The experimental realisation of a weakly interacting BEC in a magnetic trap achieved in 1995 by E. Cornell 
and C. Wiemann at Boulder and W. Ketterle at MIT in vapors of 87 Rb (a s = 5.77nm) and 23 Na (a s = 
2.75nm). This became possible by a combination of evaporative cooling and laser cooling. These systems 
are well described by Bogoliubov theory and the GP equation. 

For models of the trapped condensates as those realized in experiments, usually a harmonic trap poten- 
tial of the general form 



Vext(r) = V tI (r) = -{u 2 x x 2 + Lo 2 y y 2 + u; 2 z 2 ) (14) 



m 
~2 

is assumed. For an ideal Bose gas, the critical temperature is given as [4| 



in contrast to the critical temperature of a homogeneous BEC in Eq. (Q}. Instead of Eq. (O, the condensate 
fraction in a trapped condensate is 



no 
ntot 



if T>T C 



In rotating BECs, quantized vortices and vortex lattices have been observed, a phenomenon which is 
also known in type-II superconductors and superfluid 4 He |[T6l[T7ll . Vortices are observed by absorption 
imaging lfl8l . 
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If the condensate is in rotational equilibrium at angular velocity £1 around the z-axis, the critical angular 
velocity fi c , at which the creation of a vortex occurs, as well as the stability and dynamics of vortex 
cores and vortex lattices have, can be calculated by minimizing the free energy within the GP approach 

USES EDI221 El- 



1.3 Light scattering and structure factor 

Light scattering experiments on BECs allow the study of density fluctuations. In so-called Bragg scattering 
experiments, light scattering is studied as a stimulated process, induced by two laser beams which illumi- 
nate the atomic sample [24]. In scattering events elementary excitations are created, and the momentum 
and energy transfer is pre-determined by the angle and frequency difference between the incident beams. 

The most important quantity here is the dynamic structure factor »S(q, ui), which is proportional to the 
excitation rate per particle. Here, q = qt — q;, and q; is the wave vector of the incoming, qt is the wave 
vector of the reflected light beam, and ui is the frequency difference between the two laser beams. 

The dynamic structure factor describes a correlation between a density fluctuation at time to = and at 
time t\=t and is defined as the expectation value l25ll 

S(q, W ) = -!- / (p q (*)p+(0)) e^dt , (17) 

-''tot J 

with the density operator in momentum space, which is given as 

f>t= I ^e^ r d d r = ^a+ +q a k , (18) 

J k 

in Schrodinger representation and 

£ q (t) = e -i(*-M*)t//i^ ^{H-^t/h (19) 

in Heisenberg representation, where and a^, a£ fulfil bosonic commutation relations. Integrating over all 
frequencies u> one obtains the static structure factor 



S(q) = J S(q,w)du;, 



(20) 



which is equivalent to the line strength of the Bragg resonance. The static structure factor is then given by 
Eq. © as 

S(q) = ^- (/5 q (0)p+(0)) . (21) 

i 'tot 

In the ground state of a non-interacting condensate, the static structure factor is unity, and in the Bogoliubov 
ground state, it is given as 

h 2 a 2 

S «" " 2^k • (22> 

where _E q is the quasiparticle spectrum given in ([8]). This result has been originally derived by R. Feynman 
for the static structure factor of superfluid 4 He (26J, and will be reproduced in chapter [4] In the regime of 
long wave lengths this becomes 



where c is the sound velocity. 
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Fig. 1 Zero temperature phase diagram of the Bose-Hubbard model calculated in mean-field theory. 



1 .4 Optical lattices 

Recently, ultracold gases were superimposed by optical lattices, which are created by standing waves of 
laser fields ll27l . There are one-, two- and three-dimensional optical lattices. The lattice potential of a 
three-dimensional cubic optical lattice created of three perpendicular laser beams parallel to the coordinate 
axes, is of the general form 

Matt (r) = V x sin 2 (q x) + V y sin 2 (q y) + V z sin 2 (q z) , (24) 

where the amplitudes V x , V y , V z are proportional to the intensity of the laser field. Together with the 
harmonic trap potential given in Eq. (fT4l i the external potential of the atoms is V cx t (r ) = Vt r (r ) + Vi at t (r ) . 

A one-dimensional Bose gas, where the movement of atoms is only possible in one direction (e.g. the z- 
direction), can be created by tightly confining the particle motion in two directions (the x- and y-direction) 
to zero point oscillations. This can be done by increasing the amplitude V x and V y until tunneling of atoms 
through the lattice wells is prohibited. If V z = 0, the Bose gas is trapped in one-dimensional tubes, and if 
V z but small compared to V x and V y , a one-dimensional lattice is created where atoms can only tunnel 
between neighboring lattice-sites in the z-direction [28 1. 

The conventional model for a single-component system of bosons in an optical lattice is the Bose- 
Hubbard model. Assuming a d-dimensional simple-cubic lattice potential with q x = q y = q z = q and 
V x = V y = V Z = Vb/3, it has the form (29] [30] ED 

H B h = a+a r , +y"V r a+a r + ^ ^ a+a+a r a r , (25) 

(r,r') r r 

where r, r' denote the discrete positions of the lattice sites, a and a + are bosonic annihilation and creation 
operators and the sum of the kinetic term runs over nearest neighbor sites only. The position of site i is 
at a minimum of the lattice potential, i.e. V[ a tt(ri) = 0. 

The Bose-Hubbard model can describe a new phase, the Mott-insulator (MI). It is characterized by a 
complete loss of phase coherence between different lattice sites and an integer number of bosons at each 
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Fig. 2 A projection of the phase diagram of the Bose-Hubbard model in the vicinity of the point, where the two Mott 
lobes meet, n and J are in arbitrary energy units after the projection. 



lattice site ("lobes" in the phase diagram in Fig. Q]). The loss of phase coherence has been shown in 
experiments 11271 . The MI is favored if the on-site interaction U dominates the hopping J. 

In the hard-core boson model, which will be discussed in the following sections, each lattice site cannot 
be occupied by more than one boson. Contrary, the Bose-Hubbard model which allows multiple occupation 
to the price of the interaction energy U. The existence of BEC phase in the three-dimensional hard-core 
boson model has been proven rigorously ll32l . 

The Hamiltonian of the hard-core boson model can be written in terms of creation and annihilation 
operators a+ and a r with the usual bosonic commutation relations [a r , a J] = for different sites r ^ r'. 
They have the additional hard-core property 

a 2 r = (a+) 2 = , (26) 



which limits the occupation number at lattice site r to and 1. With those operators, the Hamiltonian is 

ma si 

<r : r'> r 

The hard-core boson model can be understood as a projection of the more general Bose-Hubbard model 
in the vicinity of those points of the phase diagram, where two adjacent Mott lobes meet (Fig. |2). This is 
similar to the picture which was applied to the tips of the Mott lobes in a recent paper by Huber et al. [35 1. 
It is based on the following idea. The number of bosons per site is fixed in the Mott state. For adjacent 
Mott lobes this means that the corresponding Mott states differ exactly by one boson per site. Now we 
consider two adjacent lobes with n and n + 1 (n > bosons per site), respectively and assume that the 
chemical potential is fixed such that the ground state is the Mott state with n particles per site. Low-energy 
excitations of this state for a grand-canonical system are states, where one or a few sites (e.g. k > 1 
sites) have n + 1 bosons, all other sites have n bosons. The k excessive bosons are relatively free to move 
from site to site on top of the n Mott state. Therefore, the physics of these excitations can be described 
approximately by the tunneling of the k excessive bosons alone. Due to the repulsion of order U, assumed 
to be not too small, it is unlikely that a site with n + 2 bosons is created. Consequently, these excessive 
bosons form a hard-core Bose gas. 
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1 .5 Outline of the following sections 

In section |2] the functional integral representation is introduced in the form as it is applied to the models 
which are reviewed. It is shown that all physical quantities can be derived from of the functional integral 
representation of the grand canonical partition function. 

In section [3] exactly solvable models are presented, namely the ideal Bose gas and a one dimensional 
hard-core Bose gas an optical lattice. Section|4]presents a summary of the results of the weakly interacting 
Bose gas on the level of Gaussian fluctuations around the mean-field solutions. It leads to the well-known 
results of Bogoliubov theory. Two approaches to the dense regime of strongly interacting bosons are 
provided in section|5]The first one will be called the paired-fermion model, and the second is based on the 
slave-boson approach. 



2 Functional integral method 

2. 1 Grand canonical partition function as functional integral 

The grand canonical partition function Z of a many-body system contains all information about the ther- 
modynamic equilibrium properties of that system Q. For given Hamiltonian H it is given as the trace of 
the density operator p: 

p = e -ftB-rf«*) . Z = Tr (p) (28) 

Here, (3 = l/(k-QT) is the inverse temperature, p is the chemical potential and the particle number operator 
is TVtot = Y^a ®aQa- It is possible to write a grand canonical partition function in terms of a functional 
integral [36] 031 • 



2.1.1 Bosonic functional integral 

Consider a bosonic many -body system given by the Hamiltonian i?(a+, a a ), where the creation and anni- 
hilation operators a+ and a a fulfil bosonic commutation relations: 



[a a , Og]_ = S a/3 ; [a a , ap\- = [a+, a^]_ = 



(29) 



The index a denotes the states |a) of an arbitrary single-particle basis, e.g. a can denote a lattice site or a 
wave vector. The grand canonical partition function is given as a functional integral over the complex field 



Z = lim I i 

M^oo 



M 

-A{<p\<t>) jj jj: 

n— 1 a 



2ni 



with the action 



£ 

M 



M ( 
n — 1 ^ a 



M 



We require for bosons the periodic boundary conditions <f> a \ 



and i 



(30) 

#(C,„+i,</w)j ■ (3D 
i = 4>* a M+v The function 



H(4>a n+ i, 4> a ,n) is obtained from the Hamiltonian H(a+, a a ) by making the replacements a+ — > <p* a n+1 
and a a — * cf> a ,n- After performing the limit M — > oo, n plays the role of a continuous imaginary time 
variable. Using r := nh/3/M we can write 

Z = J e-^^> Wt)<Ht)) , V(r(r)Hr)) := Jfa^ f[ ]J (32) 
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and 

A(0%^) = ^ W dr|^0*(r)^| : - At )^(T)+if(0*(r) ) ^(r))| . (33) 

In the following we keep M finite during the calculations and the limit M — ► oo is performed in the end. 

2.1.2 Fermionic functional integral 

In the case of a fermionic many-body Hamiltonian if (c+, c a ), the creation and annihilation operators fulfil 
the anti-commutation relations 

[c a ,cp}+ = 8 aP ; [c a ,cp] + = [c+,c^]+ = . (34) 

A functional integral of a fermionic system is given as an integral of conjugate Grassmann variables. The 
definition of a Grassmann algebra can be found in refs. 1361 [T4l [371 . Here it shall only be mentioned that 
the variables of conjugate Grassmann fields ^, tp are anti-commuting, i. e. 

and a Grassmann integral gives unity only if it is performed over a full product of all variables, and zero 
otherwise: 

i>a,ni>a,ndi) atn di> atn = 1 , (35) 



dip a . n dip a . n = / ip a .n dip a . n dip a . n = / ip a .n dip a . n dip a . n = . (36) 

Using these rules, the functional integral of the fermionic grand canonical partition function can be 
constructed by analogy with Eq. ( f30b as 



M 

Z = lim / e - A ^ TT TT d^, n d^,n ■ (37) 



7i—l a 



In the action OTI ), the complex variables 0* „, <j> a ,n have to be replaced by the Grassmann variables 
ipa,n,4>a.n, and the periodic boundary conditions have to be replaced by anti-periodic boundary condi- 
tions %jj a .i = —if> a ,M+i and -0 Q .i = -^a if+i' The same replacements can be done in the imaginary time 
functional integral defined by Eqs. (l32l and d33l , then the integration measure in d32b is replaced by 



M 

V$(t)iI>(t)) := lim TT TT d^ Q ,„d^ Q ,„ (38) 

_/V/ — » run A A A J. 



for the Grassmann fields. (For the construction of the functional integral for bosons and fermions with 
coherent states see Appendix iBli 

2.2 Correlation functions 

Physical quantities can be written in terms of expectation values. The expectation value of an arbitrary 
operator X is given by the relation 

(X) = ±Tr(xfi) (39) 
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with the density operator d28l >. The general static n-particle correlation function (CF) is defined as a product 
of n creation and n annihilation operators: 

C„(ai, ...,a n ; /?„,.. . ,/3i) := (a+ • • • a^a^ • ■ • a 0l ) . (40) 

In the functional integral representation of a bosonic system, an expectation value of some function 
f(4>*, <j>), which depends on the complex field variables, is defined as 

(/(**, 0)) = \ J f(P,4>) e~ A ^^V{r{T)<P{T)) . (41) 

Note that in a fermionic system, the complex fields have to be replaced by Grassmann fields, otherwise 
there is no difference in the formalism. To translate the static CF (l40l > to an expectation value in terms of a 
functional integral, it is necessary to introduce a dynamic n-particle CF, which depends on the imaginary 
time variable r. Therefore we introduce the imaginary time Heisenberg representation of the bosonic 
creation and annihilation operators a+ and a a : 

a a (r) = e^-^y h a a e- T ^-^^/ H 
The dynamic n-particle CF can now be defined as 

C n (am, . . . , a n T n ; PnTn+l, ■ ■ ■ ,fhT2n) '■= (^(n) ' • ' ( T n)a/3„ ( T n+1 ) ' • ' Oft (r 2 „)) . (44) 

An expectation value of the complex field variables is given as an expectation value of a time ordered 
product of the creation and annihilation operators in the Heisenberg representation |36|. The time ordering 
in the imaginary time variable is indicated by the time ordering operator T, The ordering begins with the 
largest imaginary time and ends with the smallest. The rule for a translation of an expectation value of 
a time ordered product of operators into an expectation value of a product of complex field variables is 
simply 

(d( T l) ' ' •C„( T n) < /W + l( T n+l) ' ■ ■ <Pa 2n ( T 2n)) = 

(fa+(ri) • • ■ a+ n (T n )a an+1 {T n+1 ) ■ ■ ■ a a . 2n (r 2n )) . (45) 
Introducing a time-slice e > 0, the static n-particle CF (l40b can thus be constructed by 

C n (ai, . (3 n ,...,(3x) = 
lim (ajj (r + (2n - l)e) ■ ■ ■ a+ n (r + ne)ap n (r +(n— l)e) ■ ■ ■ ap 1 (r)) = 



(42) 
(43) 



lim {<f ai (r + (2n - l)e) ■ • • 0^(r + ne)^ n (r + (n - l)e) • • • ^ (r)) (46) 

Note that this expression is independent of r. Because the imaginary time is periodic with periodicity h{3, 
it does not matter which point r is regarded as the beginning of a period, thus in particular we can assume 
t = 0. In general, it is not possible to replace the limit e — > simply by putting e — 0, because the 
limits for s > and e < are not necessarily the same. This feature reflects the fact that the creation and 
annihilation operators do not commute in the operator formalism. 

Some relevant physical quantities which can be calculated from correlation functions shall be mentioned 
here: 
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2.2. 1 Total particle number 
The total particle number is derived from the grand canonical partition function by [3 ] 

Applying Eq. d47| > to Z as it is given in Eqs. d32b and d33l , we get 



(47) 



lim 



1 1 



E 



hi) 



+op Z 

Because of the independence of the CFs of r, we have 

iVtot = lim £>*(e)&,(0)) . 

ot 

The particle occupation number in state a is 



n a = lim ((f)* (e)(/> Q (0)) . 

£^0 



C(t + £)</-« (r)dr 



-^W".«2?(^*(t)^(t)) 



(48) 



(49) 



If a denotes a position in space or a lattice site, n Q is a local particle density, if a is a momentum index, 
n a is the momentum distribution of particles. 

As has been mentioned before, it is not allowed to put the time-slice e = in general, because in the 
discrete-time definition of the action (l3TT l. the /i-dependent term is given by 



£ 

M 



M—l 



(50) 



7i— a 



and therefore occupies the off-diagonal matrix elements in the imaginary time index. It should be noted 
here, that it is also possible to construct the functional integral with the /i-dependent term being on the 
diagonal matrix elements, i. e. 



TIZZY. ^C,n^,n 



M—l 



M 



(51) 



In this case the occupation number would be n a = (0* (0)<f) a (Q)), which means that the expressions for 
the physical quantities significantly depend on the definition of the functional integral, which in some cases 
might be more convenient. However, in this chapter we will keep the off-diagonal representation given in 
(SB. 



2.2.2 Condensate density 

The condensate density of a BEC is a measure for the off-diagonal long range order of the one-particle CF. 
It has to do with the spacial range of the one-particle CF and thus a should denote a position vector (in a 
continuous system) or a lattice site (in an optical lattice). In terms of complex variables, the definition ^ 
of the condensate density in a system without confining potential is 



no := lim lim(cj)*(e)4> r '(0)) 

r r' — >oo £ — >0 



(52) 
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2.2.3 Density-density correlation function 

The density-density CF is a two-particle CF. It describes the spacial behaviour of density correlations, 
which means that here a denotes a position index as well. In terms of field operators it is defined as 

D(r - r') = (h r n r ,) = (#> r ^ r ,> , (53) 

and in terms of complex field variables it is given as 

D(r - r') = lim («(e)0r(O)#,(e)MO)) (54) 

A good physical quantity, which describes correlations of density fluctuations is the truncated density- 
density CF 

Arunc(r - r') = (h r h r >) - (n r )(n r >) . (55) 
The Fourier transform of the density-density CF is called the static structure factor 

S W = aT" £ D ( r - r')e ik ' (r - r,) • (56) 



tot 



3 Exactly solvable models 

3.1 Ideal Bose gas 

3.1.1 Hamiltonian and partition function 

In this chapter we will survey the basic results of the previously mentioned quantities for an ideal Bose gas. 
This seems to be reasonable, because it allows us to introduce the methods we will apply for an interacting 
hard-core Bose gas as well. Contrary to the interacting system, exact analytic results can be found for the 
non-interacting case of the ideal Bose gas. 

A non-interacting Bose gas in a rf-dimensional cubic lattice with nearest-neighbour hopping J and 
lattice constant a is given by the Hamiltonian 

H = J — ^— ^2 a+a r i + J ^2 a^flr (57) 

(r,r'> r 

with the dispersion relation 

J d 

£k = J — -: cos(afc„) , (58) 

where k u is the ^-th component of the rf-dimensional wave vector k. Note that the sum over nearest 
neighbors (rj, rj) means, that the index i runs over the entire lattice and the index j runs over all sites, 
which are nearest neighbours of j. This means, that each bond appears twice in the sum, once with a 
hopping process from site i to site j and vice versa. For small wave vectors k, the lattice dispersion can be 
approximated by the translation invariant counterpart 

(k =iA + ( k *), m -:=%, (59) 



where m* is the band mass. 
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We apply the discrete time action given in Eq. OTt and perform the limit M — > oo at the very end. It is 
possible to write the functional integral OTb in the form 



Z = Km / exp 



M 



E E <&,„^k. 



M 



k n=l 



k n,m=l 

where the relation between the complex fields in real space and in momentum space is 



^ k 



(60) 



(61) 



and the matrix elements of represent the structure of the discrete imaginary time variable: 



1 










-6 k 


-6k 


1 












-6k 


1 











-6 k '■ 














1 












• -6k 


1 



6k = 1 - -^(e k - n) 



(62) 



The entry in the upper right corner is necessary to realize the periodic boundary conditions. The Gaussian 
integral can be integrated out and we get 



Z= lim TT det i (k) = lim TT 



1-1- 



g(fk - m) 
M 



M 



If we now, as a final step, perform the limit M — ► oo, we get the correct form of the grand canonical 
partition function of an ideal Bose gas [36|: 



Z = H [l - e -0C«k-ri 

k 



(63) 



3.1.2 One-particle correlation function 

As already discussed in section [2721 the momentum distribution and the condensate density in a Bose gas 
can both be described by the one-particle correlation function, cf. Eqs. d48l l and ( |52l . Thus we should at 
first calculate the one-particle CF for an ideal Bose gas in general to determine those quantities. To achieve 
this we again start with the discrete time functional integral and take the limit M — > oo at the end of the 
calculations. In this sense, we define the imaginary time dependent one-particle CF in momentum space as 

C(ki,Ti;k 2 ,r 2 ) = {<& uni <t>k a ,n3) = 



lim — 

iW-xx> Z 



9k 1 ,ni9»ka,T. 2 exp 



M 



E E <t>l,n A( nl<P^ 



k n 7 m—l 

where the indices ni, are defined such that 

^K 2 -l)<ri, 2 <-r 11 ,2. 



nn 

k n=l 



(64) 



(65) 
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The Gaussian integral d64b picks out a matrix element of the inverse matrix A 1 

)n 2 ,«l ^ k l> k 2 • 



Cfc.Tijka.T*) = lim (A^)-i <5 kl ,k 2 • (66) 



M— too 

Therefore it is necessary to determine the matrix elements of A -1 . By means of the unitary transformation 
matrices 

U = 1 e^r" m f/ + = 1 e~^r" m (67) 

V M V -"2 

we can diagonalize the matrix to get 

(u(AM)-iu+) jn = **> , 

1 — Oke m n 

(i (k) )^ = p^(cr(i«)-i^ )I7 ] Jfl = x: ^ a : . 

This sum is given in the Appendix. The result is 

fi(k)\-i _ _J_ v / 4"" if n /fiON 
1 - o^ J I °k if J < n 

Performing the limit M — » oo in d65b we get 

<5 kl ,k 2 f e (^-n)(^-M)/fi if Tl > T2 



C(ki,ri;k 2 ,T 2 ) = - - e _^( e 2 k - M ) x ( efa-n-WXek-M* if n < T2 ' (69) 
Using this result and the definition (l46*l l. the one-particle CF in momentum space for an ideal Bose gas is 

d(k; k') = lim (<&(e)«M0)) = _ 1 = «k,*JVk , (70) 

where is the usual momentum distribution of an ideal Bose gas. 

In the condensed phase, where the chemical potential takes the value /i = 0, the momentum distribution 
function diverges at k = 0. In this case, the lowest momentum state k = is macroscopically occupied 
and builds the condensate. The condensate density in this case is given by 

no = w • (71) 

The normaliation with the number of lattice sites TV is necessary, because in the BEC phase the ground 
state is the only macroscopically occupied state, whereas all other occupation numbers are of the order 
of unity. The total particle density in the condensed phase is the sum of the condensate density and the 
particle density of all excited states. In the thermodynamic limit, the sum becomes an integral: 

f d 3 k 

That = tiq + J N k . . (72) 

It should be noted here, that in one and two dimensions a condensate cannot exist. The reason is, that the 
integral (|72l is divergent in these cases if [i = 0, because rtk behaves like ~kr 2 for small momenta. 

This definition of the condensate density in an ideal Bose gas is also compatible with the more general 
definition via off-diagonal long range order given in Eq. d52| >: 

lim lim<#(e)4v(0)> = lim C(r;r') = lim [ A^N k ^ r -^ 
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3.1.3 S true ture factor 

From Eqs. d54l > and d56l > the static structure factor can be obtained. The fourth-order correlation function 
can be calculated using Wick's theorem (AppendixO: 

^(«(e)W)^( £ )^(°)) = 

lim[<^( £ )0 k+q ( £ ))<^ +q (O)^(O)) + W(£)^(0))(^(0)Ue)>] = 

^k<5k,k+q-/Vk'^k'+q,k' + ^k^k,k' (Mc+q + 1) <!>k'+q,k+q • 

For q ^ 0, the first term vanishes. Thus we find the result 

^) = ^E Ar k(A r k +q + l). (73) 

k 

In the BEC by separating the ground state and excited states, we get 

5(q) = 1 + 2n 7V q + _L V 7V k iV k+q . (74) 

k#{0,-q} 

Instead of Eq. d54b one can use the more convenient definition in terms of expectation values without time 
slices 

5(q) = 1 + -!- E <^k(O)0k' + q(O)0k +q (O)0k'(O)) , (75) 

Ntot k,k> 

which leads to Eq. d74l i as well. Graphs for different temperature regimes are shown in Fig. [3] 

3 . 1 .4 Random walk expansion and world-lines 

In this section a very intuitive method of diagrammatically visualizing a grand canonical partition function 
shall be introduced for an ideal Bose gas in an optical lattice, namely the random walk expansion Il38ll391 . 
We will perform the same expansion in the following chapters for a system of hard-core bosons, in order 
to demonstrate the effect of the hard-core condition. 

The grand canonical partition function of an ideal Bose gas in a d-dimensional cubic lattice is given by 
the functional integral Eq. d60l i. But here we use the real-space representation. The time structure of the 
matrix A is the same as in Eq. d62b , but instead of the dispersion relation e k we use the hopping matrix 

- J —J /Id if r, r' nearest neighbours 

Jrr ' \ otherwise ' ' ( 6) 

which establishes the spacial structure of A, and make use of 

e rr / := J rr / + J8 rT > . (77) 
Thus we can write 



(78) 



where the term d n iS m M accounts for the upper right matrix element in (f62t which arises from the period- 
icity in imaginary time. 
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S(q) 



T = 

T = 0.1 T c 

T = 0.5 T c 

T = T c 

T = 1.25 T r 



\ 



0.0 



0.5 



1.0 



1.5 



2.0 



Fig. 3 Static structure factor of an ideal Bose gas of free particles. At T — 0, S is constantly unity and has a 5-peak 
at q — 0. At < T < T c it diverges, and at T > T c it reaches a constant near q = 0. All cases are characterised by 
the relation lim^oo S(q) = 1. 



The idea of the random walk expansion is to expand the off-diagonal part of the exponential in the 
functional integral expression in terms of the field variables: 



exp 



M 



<Pr,n^-rr':nm(f>r 

r,r' n,m=l 



exp 



A I 



r n=l 



y — 

. L rr / .n- 



-I 0„',n>0} 



r,r',n ^ 



M 



(e rr / - /U^rr') 0r',n-l 



The abbreviation M rr ' has been introduced for convenience. The functional integral can be solved by using 
the identities 



II (<n<n-l) irr ' n =n[(<") m '' n (^) <n ] > 
r,r',n r,n 

where m r , n := ^ i rr /, n and mj. >n := ^ k'r,n+i 



(79) 



and 



- — III'. <J mrn ' 



2m 



(80) 



This results in the following form of the grand canonical partition function as a sum over all indices l rr >, n '- 



{irr',„>0} r >™ <■ 



■)n 



fcrr')'"'- 
l Yr' n ■ 



(81) 
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Fig. 4 Random walk expansion of an ideal Bose gas; world-line diagram. 



Note that it is necessary to define (w rr ')° = 1 here, even for the vanishing matrix elements of u. 

One possible interpretation of this expression is as follows: Each term of the sum can be represented 
by a diagram, where a particle propagation from site r at imaginary time r to site r' at time r + h/3 /M is 
indicated by an arrow. So each particle is characterised by a "world-line" showing its movement through 
the lattice in imaginary time. The contribution of a certain diagram is defined by the following properties: 

• The number of particles (arrows) propagated from site r' at time (n — I)hf3/M to site r at time 
nh/3/M is given by l rr >.n- In the case of nearest neighbour hopping, particle propagation in one time 
step h(3/M is only possible between neighbouring sites, or the particle stays at the same site. 

• The number of particles (arrows) which are propagated to site r at time nh(3 /M from the previous 
time step is m ri „. 

• The number of particles (arrows) propagating from site r at time nh(3 /M to the next time step is m' T n . 

• Particle conservation is assured by the (5-function in Eq. (ISTI l, such that m r „ = m' r n is equal to the 
number of particles at site r and time nh/3 /M, 

• There is a periodicity in imaginary time: Time t = h/3 is equivalent to time r = 0, so the diagrams 
have to be periodic in time. 

Note that in the ideal Bose gas m r .„ > 1 is possible, i.e. more than one particle can occupy the same lattice 
site at the same time. This will be excluded to establish the hard-core interaction in a Bose gas. 

3.2 Hard-core bosons in ID 
3 .2. 1 General remarks 

The main feature of the one-dimensional hard-core Bose gas is, that the particles cannot interchange their 
position. An interesting consequence of this property is the equivalence to an ideal non-interacting one- 
dimensional Fermi gas. However, it is important to mention, that this equivalence does not hold for all 
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physical quantities in momentum space, namely for those which are given by one-particle correlation 
functions like the momentum distribution 140, 41, 42 43 44 1 . It is possible to calculate the momentum 
distribution by means of a Jordan-Wigner transformation (see e. g. refs. [36, 45]). This approach has been 
used in a couple of works 0281 1461 l47l . However, this problem will not be addressed here. On the other 
hand, quantities given by two-particle correlation functions like the density-density correlation function 
and the dynamic structure factor are the same for hard-core bosons and for ideal fermions. 

The zero temperature phase diagram of a hard-core Bose gas in a one-dimensional optical lattice shows 
three phases l48l : An empty phase (EP), an incommensurate phase (ICP) with a particle number per lattice 
site of < n to t < 1, and a Mott insulator (MI) with n to t — 1. Here we will especially be interested 
in the phase transition between the ICP and the MI phase for zero and non-zero temperatures. Again, the 
quantity we chose for investigating this transition is the static structure factor. It has also been considered in 
other works about one-dimensional Bose gases, in the weakly interacting regime as well as in the strongly 
interacting regime l|49l l50l IBT1 l52l . 

As has been demonstrated for the ideal Bose gas, a random walk expansion leads to a world-line picture. 
To make the mapping to a system of ideal fermions possible, it has to be assured that world lines cannot 
intersect each other. So instead of constructing the functional integral by starting from the Hamiltonian, 
we choose a different way and construct it by starting out from the random-walk picture directly. 

When the random walk expansion for a system of ideal spinless fermions is performed, one obtains a 
sum which is analogous to the sum in Eq. ( T8TT > with two important differences: Because of the nilpotent 
property of the Grassmann variables, the fermionic analog to Eq. ( f80b reads 

$"V m ' e"^ = S mm ,(5 m , + 6 mA ) . (82) 

This means that all terms, where the particle number m r n or m' x n is larger than 1 at lattice site r, do 
not contribute. This reflects the Pauli principle or in the case of hard-core bosons, the hard-core property. 
The second is that the Grassmann variable analog to Eq. (|79| > gets an additional sign because of the 
anti-commutation property. To avoid this problem it is possible to construct a world-line model where 
world-lines do not intersect. For this purpose we adopt an approach to the statistics of directed polymers 
in two dimensions ll53l . 



3.2.2 Particle density and phase diagram 
It has been shown that the grand canonical partition function is given by the functional integral 

M 2 



Z = lim / exp 

M->oo 



with the 2 x 2 matrix 



EE E^J^ff^ 

k n=l j,j'=l 



n,3 



JJ dlpk,n,j dlpk,n,j (83) 
k,n,j 



G- 1 (k)=[ R , + ., ^ M 2 lf, n 1 + ' . (84) 

This integral can be performed and it yields 

f g X-2A/A/- 

Z= lim 1-iiu) detG" 1 , (85) 

M-»oo \ M J 

where JV is the number of lattice sites. The one-particle correlation function of the fermions at equal times 
can be calculated as 



1 M 

C(k)= lim — V Gn(k) nm (86) 

M— *oo IVl A — ' 



n,m— 1 
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"tot 



0.5 




Fig. 6 Phase diagram of the one-dimensional hard-core Bose gas at zero temperature with an empty phase (EP), an 
incommensurate phase (ICP), and a Mott insulator (MI). 



This sum is performed in Appendix lA.4l After performing the limit M — > oo, the result is 




As was mentioned before, the one-particle correlation function does not lead to the momentum distri- 
bution. However, the total particle density of the bosons is given by taking the fermionic one-particle 
correlation function in real space 

C(r,r')= f C(k)e ik ^ r - r ">^ (88) 
Jo 2vr 



21 



at r = r'. This can be shown by applying the expression ( f4Tb of the total particle number to the partition 
function d83l >. We need an additional factor of 1/2 because of our special construction: 



Id Id 
JVtot = — —\ogZ= lim — — 



-2MAHog ( 1 - -j^/x ) - log det G 



1 P 

r,n,j 

So, because of (Vv,n,i^V,n,i) = ffr,n,2^n2) = C(r, r), we find the result 

ntat = ^ = 1-C(r,r) (89) 

for the total particle density. Note that the time slice e, which was necessary for the definition of the total 
particle density for weakly interacting bosons (see Eq. d48l>). is absent here, because of the construction of 
the Green's matrix. The zero temperature result is 



lim n to t 



if/x<-J 

1 - i arccos (f) if —J < fi < J . (90) 
I if fi < J 



Graphs for zero temperature and finite temperature are plotted in Fig. [5] Both graphs are symmetric to the 
point fi/J = 0, ntot = 1/2. This reflects the particle hole symmetry of the system: Because of the Pauli 
principle a given configuration of the system is symmetric to the configuration, in which each occupied 
site is empty and vice versa. Further one can see that the system is empty (n to t = 0) if fi/J < — 1, and it 
is a Mott-insulator (n to t = 1) if m/J > 1. The phase transitions between the EP and the incommensurate 
phase with < n to t < 1. and between the ICP and the MI, are characterised by a diverging slope of the 
curve at the transition points. At non-zero temperatures the sharp phase transition is smeared out. The zero 
temperature phase diagram is depicted schematically in Fig. [6] 

3.2.3 Density correlations and static structure factor 
We define the truncated density-density CF of the hard-core Bose gas as 

D(r - r') = (^r,n,lVV,«,lVV',n,lVV,n,l) ~ (VV,n,lVV,n,l) (VV',™,lVV,n,l) ■ (91) 

V v ' 

= "tot 

Using Wick's theorem for Grassmann variables as given in Appendix|C] we find 

(Vv.n.iVv^.iVV' ,n,iVV' ,71,1 ) = "tot - C(r,r')C(r' , r) , 

leading to the result 

D(r-r') = -C(r,r')C(r',r) . (92) 

The static structure factor is related to the density-density CF by means of a Fourier transformation which 
is shifted by unity, and a normalisation. We use the definition lt48l l50l 

S(q) = 1 + n( ,\ ■ (93) 

It is the analog to the definition of the static structure factor of an ideal Bose gas ( f75l l, where the term 1 
appears when the time slice is canceled in the expectation value of the complex fields. Expressed in terms 
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of the one-particle CF in momentum space C(k) by applying the Fourier transformation in Eq. ( 1881 ), the 
above expression reads 

We want to investigate the static structure factor at zero temperature in the ICP near the phase transitions 
to the EP and the MI. Because of the particle-hole symmetry discussed in the previous section, both tran- 
sitions should be symmetric with respect to the physics of light scattering. Let us first discuss the region 
/i > 0. Defining the characteristic wave vector k* we find the result 

r ^ if(z<2fc* 

S(q) = \ 1 if 2k* < q < 2tt - 2k* . (95) 
I ifq>2n-2k* 

In order to keep the particle hole symmetry for the static structure factor, in the region [i < we make 
the substitution C(k) — ► 1 — C(k) in the expression d94l ), and find the same result as in Eq. ( l95l l. The 
expression for the density-density CF D(r — r') near both phase transitions we get from the Eqs. d88l and 
d92l ). At zero temperature it is 

v 7 V 2it(r-r') / 
The characteristic wave vector can be written in terms of the total particle density fl89l >: 

* f 27rntot if "-tot < 1/2 q? » 

" \ 27r(l-n tot ) ifn to t > 1/2 ' 1 ; 

Near the phase transitions where 5 := \[i — /Lt c |/J <C 1, we have fj, = (1 — <5)J at the ICP-MI phase 
transition, and = — (1 — S) J at the ICP-EP transition. Here, we can approximate 



k* « V8J . (98) 

For a homogeneous impenetrable Bose gas the role of k* is played by the Fermi wave vector fcp = 7rntot 
1 50]. In our result ((97}, fc* depends linearly on the density as well as in the region n to t < 1/2, but the 
discontinuous slope of the function fc*(ntot) at the point ntot = 1/2 is a consequence of the optical lattice 
potential. The relation (l23l l allows us to identify the excitation spectrum 



frk* 

e(q) = hcq + 0( q 2 ) , C= . (99) 

m 

which is linear for small values of q, where c is the sound velocity. The density-density CF and the static 
structure factor near the ICP-MI phase transition are plotted in Fig. [7] 

The density-density CF shows characteristic oscillations with length A = ir/k*. This length scale 
diverges at the ICP-EP and ICP-MI phase transition with l/n to t and 1/(1 — Utot)> respectively. Thus it 
can be used as a measure for the distance of the system to one of the two phase transitions. In the EP and 
the MI phase, the density-density CF vanishes because of the absence of particle number fluctuations, and 
the static structure factor saturates to S(q) = 1. 

3.2.4 External trap potential 

In the previous sections a system in a translational invariant lattice was considered. Calculations have also 
been made for a one-dimensional Bose gas in a harmonic trap potential [48 1 



V(r) = -<(ar) 2 , (100) 
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Fig. 7 Truncated density-density correlation function D(r — r') and static structure factor S(q) in the vicinity of 
the ICP-MI phase transition. The transition point is at fi c = J. For the ICP-EP phase transition, the situation is 
symmetrical. 




Fig. 8 Local particle density for system in harmonic trap potential (/i = 0.7, ma 2 ui^ /2 — 3 x 10 5 ) with varying 
tunneling rate J. A Mott plateau appears in the center of the trap (r=0) as J is decreased below a critical value 
J P w 0.70. (Fig. taken from ref. (48l .) 



where again a is the lattice constant, and u\ 10 is the harmonic oscillator frequency of the trap. The numerical 
result for the local particle density at zero temperature is plotted in Fig. [8] where the formation of a Mott 
plateau can be seen below a critical value Jp. A similar behavior was found for the one-dimensional 
Bose-Hubbard model with a harmonic trapping potential |54| . 

The properties of the density-density CF and the static structure factor are qualitatively the same as in 
the translational invariant case. D(r) vanishes when Jp is reached, owing to the fact that there are no 
density fluctuations within the plateau. The characteristic length scales become larger as the Mott plateau 
is reached. 
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4 Weakly interacting bosons: Bogoliubov theory 

Before discussing an interacting Bose gas in an optical lattice, we begin with the derivation of the Bo- 
goliubov approximation for a dilute homogeneous Bose gas. Although the Bogoliubov theory can also be 
applied for bosons in a lattice potential, a Mott-insulating phase is not found within this approximation 
|55|. Many aspects of the physics discussed in this chapter show up in the hard-core Bose gases in optical 
lattices as well. 



4. 1 Derivation from saddle point approximation 

It might be interesting to derive the results of Bogoliubov theory from the functional integral point of 
view. The method which will be used here and in the following chapters is the saddle point approximation 
(or: stationary phase approximation, Gaussian approximation) ll36l l56l l57l . It allows to find a mean- 
field solution plus fluctuations around the mean-field result. The mean-field solution is connected to the 
condensate order parameter, while the fluctuations contain the information about the quasiparticles and 
their spectrum. The saddle-point approximation is good as long as these fluctuations are small. 

The main idea of a saddle point approximation is to expand the action of the system around its minimum 
up to second order in the field variables. This leads to a Gaussian integral which can be performed. The 
action of a bosonic system is given in Eq. d33i l. where in this case the index a shall denote the position 
vector r. Together with the Hamiltonian (0 of the interacting Bose gas we have 



A{4>*,4>) 



no 



dr / d 3 r{ <p*(r,r) 



8 



2m 



V 2 + V^t(r) 



(r,r) 



-f I0(r,r)| 4 



(101) 



By minimising A with respect to the complex fields we get a mean-field equation for the condensate order 
parameter <&o( r , T )'- 



2m 



V 2 + y cxt (r) + g |* (r, r)| 2 *o(r, r) 



d_ 



A» $o(r,r) 



(102) 



After performing the analytic continuation — > and omitting the chemical potential term, this 

is identical to the time-dependent Gross-Pitaevskii equation (fT~2t . We recall that the invariance of the 
mean-field solution under a gauge transformation (|9]l with global phase a reflects the broken global U(l) 
symmetry of the BEC phase. 

To find the results from the previous sections in this chapter, we assume a homogeneous system, i.e. 
T4 x t(r) = in the action dlOU . Further we assume that the mean-field solution is constant in space and 
imaginary time: $o(r, r) = $o- In this case, the solution of Eq. ( 1102b is 



l*o| 



A* 

n = - 
9 



We now write the complex field as the sum of the mean-field solution plus fluctuations 

<£(r, t) = $ + S(f>(r, t) , <f>* (r, r) = <S>* + 6</>* (r, r) , 



(103) 



(104) 



where the complex field of fluctuations 5<fi is considered to be small, such that those terms in the action 
which are of higher than second order in the fluctuations, can be neglected. We split the quasiparticle field 
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into its real and imaginary part and write 8(f>(r, r) = 8(j>' + i6<f>", 8<ff (r, r) = 8(f)' — \8(f>" . The expansion 
yields 

where we have already eliminated the condensate order parameter by Eq. 1103) . and the zeroth-order part 
of the action is 

A = /3^(-^|$ | 2 + ||$ | 4 ) =~^f- (106) 

Because Aq does not depend on the field fluctuations, and the second term is of second order in 8<f) and 
8(f)*, the functional integral for the grand canonical partition function 

e- A{s *' M " ) V{8^>'{v, T)S<p"(r, r)) (107) 

can be performed because it is Gaussian. We Fourier transform the field of fluctuations with respect to the 
spacial coordinate like 

^'(r,r) = -^=2^(r)cos(k.r) (108) 

k 

ty"(r,r) = -^=5^50£(r)cos(k.r), (109) 

k 

with the constraints 8$^ = &£'_ k and 8<p'^ = <50" k and thus get 



with the free -particle dispersion relation = fi, 2 k 2 /2m. It is further possible to perform a Fourier trans- 
formation with respect to the imaginary time coordinate as well, namely 

8&(t) = l^^cosKr) (111) 



(112) 



with the Matsubara frequencies for bosons u> n = 2im/h(3 and the constraints Scj)^^ — 8<fi' w _ u and 



u> = ^k ■ This leads to the form 



*-* + £(Sfc)-^>(Sfc)- 

k,n ' ' 

and allows to identify the quasiparticle Green's function (a 2 x 2 matrix in this case) 

o- l ^)-(Z. <114) 

The excitation energies of the quasiparticles are given by the poles of the quasiparticle Green's function 
ll36ll . which are found by solving the equation 

deta _1 ( k ! i? ^n)=0. (115) 
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After performing the analytic continuation \huj n — ► we get 



#k = V e k( 2 ^ + e k ) , 
which is identical to the Bogoliubov spectrum, if the relation n = fi/gis inserted. 



(116) 



4.2 Partition function and condensate depletion 

To find the correct expression for the grand canonical partition function as well as for the correlation 
functions, we have to perform the same steps as in section 13.1.11 namely to start with the discrete-time 
functional integral and sending the number of time steps M to infinity at the end. By analogy with Eq. 
( f60b . the discrete-time version of Eq. (II 10b is 



A/ 



discrete 



k n,m— 1 



In 



k) 

)irn 



k.m 



(117) 



where A^X has the M x M structure 



i(k) 



B 


-k 








-k 


-k 


B 


-k o 









-k 


B 











-k '•■ 


-k 












B 


-k 


-K 






-k 


B 



in the imaginary time variables n and to, and each matrix entry is by itself a 2 x 2 matrix: 





Sk= K i_ ^ (ek+/x) 



1 i 
-i 1 



B = 







1 ~ m» 



(118) 



(119) 



The matrix can be diagonalized by using the same unitary transformation (1671 1. which was applied for the 
ideal Bose gas. This yields 

(UA^U+) kn = 









1-^ 



+ AO 



cos 
— sin 



s(|n) sm(fn) 



Using the product given in Appendix lA.61 the determinant of the matrix can be found as 



<i,<t + k, = (i-A (ek + M ) 



(120) 



l + 4^k( £k + 2 A1 ) + 0(^ 



71/ 



/9 



A/ 



7U 



14^(4 + 2^+0^ 



71/ 



A/ -1 



M 



(121) 



Thus we obtain the grand canonical partition function of the Bogoliubov Hamiltonian (after omitting a 
constant factor): 



Z = e 



-An 



lim 



M- 



J] [deti^ 



= exp 



k#0 



/ /3V> 2 



k#0 



cosh(/3£; k ) - 1] 5 . (122) 
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The distribution function of the particles outside of the condensate is given as 

(n k ) = <^(0)^ k (0)> = <^ k (0) 2 ) + WkW) = lim \ (([A^U + ([A^ 1 ), 

(123) 

with the 11- and the 22-component of the matrix with respect to the 2x2 structure. After inversion of the 
matrix ( 1120b and the back transformation, the matrix elements can be found and after performing the limit 

M — > oo we get 

(^) = -l+ e -^coth(^ Ek ) . (124) 



2 2E k \2 
The quantity 

/d 3 A: 

is called condensate depletion. Contrary to the ideal Bose gas it is non-zero at zero temperature. 



4.3 Static structure factor 

The static structure factor is given by the fourth-order expectation value d75l l. which we used for the ideal 
gas before. We replace 0o by the order parameter $o and for non-zero momenta we replace 0k — > <5</>k- 
After splitting the fluctuations into real and imaginary part and applying Wick's theorem for real variables, 
we get a similar result as in Eq. d74l ). The difference to the ideal Bose gas is, that the anomalous expectation 
values (0k^~k) an d (^k^-k) a l so gi ve a contribution here (for simplicity we have dropped the time 
variable). The contribution of the anomalous expectation values after splitting it into its real and imaginary 
part is 

such that the static structure factor is given as 



%) = i + 2^w + |i((^ q ) + (^ q )) + Y, <^>K +q ) 

Ntnt Ntot k^R- q} 



l + 4-^L(W q ) 2 )+ J2 NXW- (126) 



tot. 



k#{0,-q} 

After performing the limit M — ► oo we find 



((<^) 2 > = lim V (q) ]n =-- + -— cotbY^K 
u 7 m^oo2 l Jn 4 4£7 q V2 q 



(127) 



If we neglect the last term in Eq. d 1 26b which is quadratic in the momentum distribution, this expression 
reduces to 

S(q) = -g- coth (^E^ . (128) 
S(q) - e q /^q- 
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To determine the type of the decay of the density-density correlations for large distances (i.e. exponen- 
tially or algebraically) at zero temperature in d dimensions in the BEC phase, we Fourier transform the 
static structure factor for small wave vectors, because they are relevant for large distances r: 

£>(r)~ [ S(q)e iqr d d g ~ / = _ e iq r d d q ~ f , ^ e iq r d d q . (129) 

W J V J J v/2( M + J)q 2 + q 4 J v/2(// + J) 

This expression shows an algebraic decay. In d = 1 the decay is proportional to 1/r 2 (in agreement with 
the result d96*l > of the one-dimensional system), in d = 2 it decays like 1/r 3 , and in d = 3 like 1/r 4 (see 
Appendix|D]i. In the empty phase, all CFs vanish completely at zero temperature. Thus, the static structure 
factor is constantly unity. 



5 Strongly interacting bosons in the dense regime 

5 . 1 Paired-fermion model 

5.1.1 Bosonic molecules of spin- 1/2 fermions 

We now introduce a model of hard-core bosons which are constructed by molecules consisting of pairs of 
spin- 1/2 fermions, as an alternative to the hard-core boson model. In order to distinguish it from the latter 
this model will be referred to as "paired-fermion model". 

A general model which was introduced to study the dissociation of bosonic molecules into pairs of 
fermionic atoms in an optical lattice was proposed in ref. l58l . It is given by the Hamiltonian 

H - fiNtat = ™ Yl B r,a & r,a ~ ^ ^ ^^'t^l^'l "M^E ' (130) 

(r,r')<T=T4 <r,r') r cr=tl 

The index a =f, J, denotes the spin. The first term describes tunneling of individual fermions with rate i 
and the second term tunneling of local fermion pairs. Similar Hamiltonians were proposed in a couple of 
works for homogeneous systems, in order to study the BEC-BCS crossover [59, 60, 61 1. In contrast to the 
lattice-Hamiltonian ( 1130b they do not exhibit a Mott insulating phase. 

Because the main interest here shall be the model of hard-core bosons, we consider the case t = in the 
following, i.e. we exclude the existence of dissociated fermionic atoms. Further we will write the index 
a = 1, 2 as superscript instead of the spin indices f, J,. We write the grand canonical partition function 
of the system in terms of a fermionic functional integral of a field of conjugate Grassmann variables as 
defined in Eq. (|37| > with the action 

n— 1 L r,cr r.a 

+ Jj J^Jrr'^n+l^l'.n^ln+l^r^n \ > ( 131 ) 
r,r' J 

with anti-periodic boundary conditions in time. Here, we have replaced ^ — > fi/2 due to the fact that 
the chemical potential is associated with the number of paired fermions (i.e. to the bosonic molecules), 
hence the factor 1/2 in front of the term which contains fj,, while in Eq. ( 1130b . iVtot is the particle number 
operator of single fermions. 

In the world-line picture, the paired-fermion model given by Af erm is represented by pairs of fermions 
with opposite spin 1 and 2 whose world-lines always stay together while they tunnel through the lattice. 
Tunneling of unpaired fermions does not exist. The world-lines of two fermions of species 1 and 2 always 
stick together while tunneling through the lattice. 
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5.1.2 Hubbard-Stratonovich decoupling 

The idea of a Hubbard-Stratonovich transformation is to decouple a quartic term of a many-body system 
by writing it in terms of a Gaussian integral [62|. The original field variables are then only of second order 
and can be integrated out such that the system is represented only by the field variables of the Gaussian 
integral. 

We perform a Hubbard-Stratonovich transformation on the system of paired fermions |58 1 given by Eq. 
( I13U . Only the term which describes hopping of fermion pairs is quartic, so we will decouple it. Contrary 
to the case of the hard-core boson model, it is not necessary here to decouple the entire off-diagonal term, 
because the term describing the discrete-time derivative and the term containing the chemical potential are 
already of second order. For the matrix with fermionic boundary conditions we write 

-form P * 



M 



and insert the identity 



M 

const. X exp \ -JjYl n+lV& ,n^r,n+l^r 

r.r' n,m—l 



I ex P l ~ Jj ^U^rr™m) Vr'.m - ^ ^ X* r ,nXv.n 

" r.r' n,m v * ■ 



,n T X 



•.«)] } n 



(27Ti) 2 

(133) 



The parameter s cares for the convergence of the integral of the complex field ip. For w^r™ m we nave the 
eigenvalues 

v^ = e-^( n -i)^e k + s, (134) 

therefore one has to choose s large enough such that all eigenvalues are non-negative, but besides this con- 
dition the choice of s is free. We integrate out the Grassmann field in the functional integral representation 
of the partition function, like we did in the previous section: 

Zf CTm = / exp [-Af eim ((p ,tp,x ,X)\[[ 7^772 (135) 

r,n ^ ' 

with the action 

A{ erm (ip*,ip,X*,X) = X!X!^™^rr™m)~Vr',m + ^X!^r,nXr,n-X! l0 S det ^ ' ( 136 ) 
r,r' n,m r,n r 



where we have introduced the matrix 
r* -1 — a 



i^r.n + Xr,n 1 

1 -(i<n+X?,„) 



(§n,m+l ~ S nl S mM ) _V + 2M . (137) 
V 1 2A/ U / 
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5.1.3 Saddle-point expansion 
Under the assumption 

ft,n = </>0 <Pr,n = ¥>0 Xr,n = X* Xr,n = Xo 



(138) 



that the mean-field solution is constant in space and time, we can Fourier transform the matrix G r * = G n 1 
in Eq. ( 1137b with respect to the discrete-time index: 



g; 



l-e 



i^o + Xo 

11 2M J 



l-e M 



271/ 



By the use of the identity £ r , m (t)^ m ) _1 = (s - /3J/M)- 1 we have: 



<Po<Po . 1, 



2 ferm 

^0 _ _ 

ATM a + % ^ -5 



XoXo 



(139) 



1 M 



M 



-(i< /3o + Xo)(i^+xS)-l-e- 2 ^("-^ (l- I +2e-^( n -i) 



From the saddle point conditions 



dA fc 



dA 



ferm 



= o 



dA 



form 



dA 



ferm 



we find the mean-field equations 



= 



*° =-iG 



SI 

M 



= G 



(140) 



(141) 



(142) 



where G is calculated in Appendix [A3] and the result is 

JlPo/s 



G = 



: tanh 



J\<Po\ 



(143) 



We find a trivial solution with (fo = (pfi = xo = Xo = an( l a non-trivial solution with broken U(l) 
symmetry. For the mean-field action we find (after integrating G with respect to i(fo + xo)' 



A orm = Af 



(144) 



The complex fields <p and x are expected to fluctuate about the SP solution due to thermal and quantum 
effects. If we keep our expressions only to the first order of r = (3/M, making use of the notation 

d T = (S n ,m+i — Srum)/T and denoting A = i(j> + x and A = i0* + x*> then 



G- 1 = Go 1 



SA 





-SA 



(145) 



31 



where 



A T(d T -fi) 

r{d T + n) -A 



Applying the Taylor expansion ln(l + x) = x — x 2 /2 
logdetG -1 = trlnG -1 = trln 



Go 1 



SA 




we get 


-SA 



trlnG n 1 tr 

2 



G 



SA 
-SA 



Calculating the trace in p = {q, ui} representation we get 



-SA 



fcrm 



(146) 



(147) 



(148) 



D [Sip] exp 
where SA ie ™ is given by 

k n,m 

Here, Q represents the Green's function of quasiparticle fluctuations (Appendix[E]i. 



5 . 1 .4 Results for the paired-fermion model 

It turns out that even on the mean-field level, the paired-fermion model shows some interesting physical 
results. The condensate density we get via the definition ( f52l and the mean-field approximation that the 
CF factorizes for large distances: 



n = lim {lpl,n+l*Pr,n+l*Pr>,n*Pl>,n) = (^,n+lV'r,n+l>(V'r',nV , r',n> 



(149) 



r r — >oo 



Further, the CFs which are of second order in the Grassmann field, are given by the diagonal elements of 
the matrix G whose inverse is given in Eq. J 1 37b . These diagonal elements are equal to G/2 from Eq. 
ED: 



G 



(150) 



G 2 

2 => n ° = T- 

Thus, from the Eqs. d 143b and dl50t , together with the M — > oo limit of Eq. (11421 ). one finds a self- 
consistent equation for the condensate density: 



J = \J [i 2 + 4 J 2 n coth 



^» 2 +4J 2 n 



(151) 



The total particle density we get from the mean-field action ( 1 144-t is 



ferm 



"-tot 



1 dA% 



r 1 + ^ 



1 1 



tanh 



V" 2 + 4 J 2 



n 



I +h„il 1 ( ^ 



in the condensed phase (no > 0) 

in the non-condensed phase (n = 0). 



(152) 
(153) 



32 



Ch. MoseleJ 1 , O. Fialko 1 , and K. Ziegler 1 : Interacting bosons in an optical lattice 



1.0 



0.5 



0.0 



-0.5 



-1.0 



(a) 




MI 




^ / 

/ 

^ t 




\^ \ 
\. \ 
\. \ 
\- \ 


BEC (n + 0) 


k B T = 

k B T = 0.1 

k B T = 0.2 





0.0 0.2 0.4 0.6 0.8 1.0 

J 




-1.0 -0.5 0.0 0.5 1.0 

p/J 



Fig. 9 (a) Phase diagram with phase boundaries between the BEC and the non-condensed phase for different temper- 
atures. For &bT 7^ there is only one phase boundary between a BEC and a non-condensed phase. The energy unit is 
arbitrary because of a simple scaling behaviour, (b) Critical temperature of BEC formation. 
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Fig. 10 Total particle density and condensate density for zero temperature (thick lines, given by Eqs. fi54l and (155\ ) 
and for non-zero temperature (thin lines) plotted against chemical potential. 



It might be interesting to mention that all these mean-field results do not depend on the parameter s which 
was introduced in the Hubbard-Stratonovich transformation for the convergence of the Gaussian integral. 

The phase boundary between the BEC and the non-condensed phase we get from Eq. (115 It - The 
resulting phase diagram is depicted in Fig. [9] We see in picture (a) that for T > the phase diagram is 
separated into two parts, a BEC phase and a non-condensed phase. But at T = there are three phases: A 
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BEC, an empty phase (ntot = 0) for [i < — J, and a Mott-insulator (ntot = 1) for u > J. A density profile 
of n t ot and no is plotted in Fig. [I0]for different temperatures. At zero temperature the sharp transitions 
between the empty phase and the BEC, and the BEC and the MI, can be seen in the plot of the total particle 
density. The zero temperature result is 

n = { U 1 -^) «-J<»<J , (154 ) 
[ else 

( if/i<-J 
ntot = I 3(1-9) if-J<M<J ■ (155) 
[l if J < fj, 

If the temperature increases, the sharp transitions are smeared out. 

Calculations for the quasiparticle spectrum by finding the poles of the Green's matrix Q of the Gaussian 
fluctuations have been made for the zero temperature phase diagram ll63l . The zero temperature result in 
the empty phase and in the MI phase is 

E k = e k + H - J , (156) 

with the gap A = — J, and in the BEC phase it is 



■'''-<$)>(§)'* 



(157) 



In the dilute regime, i.e. if u = — J + A/i, with A/i <C J, this can be approximated by 



E k = ^k(2(/i + J)+e k ) . (158) 

Using the Green's function of quasiparticle fluctuations (see Appendix 10), we can calculate the effect of 
quantum fluctuations on the condensate density: 

"0 = 7 ( l - 4 )+ Sno, (159) 



4 V J 2 

where the correction to the mean-field result is 

(J 2 ~tV f d d k Blg k | (J 2 - /i 2 ) f d d k B ^ i.P -r'r 1 ,/■'*■ ul 



J 3 J (2ir) d £ k 4J 3 J (2ir) d 4J 3 J {2n) d E k 

3( J 2 - /i 2 ) 2 /* 2 f d d k B 2 <? k 



4J 5 J (2ir) d £ k ' 



(160) 



where £? k = — Y^=i cos kj> 9k.= 1 — i> k . It should be notices that this correction vanishes at the critical 

d J 

point. 

It might be interesting to mention that in the zero temperature limit near the phase transition to the 
empty phase where /i = — J + A/i with A/i <C J, i.e. in the dilute regime, it is possible to approximate 

no = ^ + G(A M 2 ) = n tot +0(Au 2 ) . (161) 



This agrees with the Gross-Pitaevskii result (1103b . if the term of order A/i 2 is neglected, and the identifi- 
cation g=2J has been made. 
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Fig. 11 Condensate density. The solid, dotted and dashed lines show the mean-field result at T — 0, the influence of 
quantum fluctuations at T = to the mean-field result and the mean-field result at T = 0.2, respectively. 

The main correction due to the thermal fluctuations are already included in our mean-field theory, where 
the condensed density is given by 



n 



4J 2 



(162) 



and \ tpo\ 2 can be determined from Eqs. (I142H1431 ). 

The effect of quantum fluctuations and thermal fluctuations is depicted in Fig. [TT| We see that both of 
them lead to a depletion of the condensate, but the quantum depletion alone does not change the transition 
points. 

The static structure factor for small wave vector q and for small temperature T in the BEC phase reads 



(J 2 - m 2 ) Jg q 



J 2 r 



E a 



coth 



f3Ea 



(163) 



where n is a total density of particles. 

In the dilute regime, i.e. close to the empty phase, when n ~ (J + fx)/J and J 



[i w 2 J we obtain 



(164) 



which is in agreement with the well-known result for the weakly interacting Bose gas (cf. section[4] In the 
dense regime, i.e. close to the Mott phase when n ~ 1, the static structure factor vanishes. 

In conclusion, we can say that the paired-fermion model has three phases at zero temperature, an empty 
phase, a MI, and a BEC, even on the mean-field level. However, at non-zero temperatures a new phase 
emerges from the MI phase and the empty phase, that is controlled by thermal fluctuations. 
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5.2 Slave-boson model 

5.2. 1 Hamiltonian and functional integral 

In this chapter it shall be shown that a slave-boson approach can be applied to describe a system of hard- 
core bosons. The slave-boson representation was originally developed for fermion systems, e.g. the Hub- 
bard model [64, 65]. It allows to account for many aspects of strong correlations even on the mean-field 
level. The slave-boson approach to hard-core bosons that will be presented here, has been developed in 
refs. [66 67, 68 69|. It is an alternative to the paired-fermion model which was discussed in the previous 
chapter. 

Again, the starting point is the Hamiltonian d27T >- We introduce bosonic creation and annihilation op- 
erators of empty (e+, e r ) and occupied b r ) sites which act on a fictitious "vacuum". To transfer the 
Hamiltonian to the extended Fock space, we replace the hard-core boson operators by 

— > b^e r ; a r — > e+6 r . (165) 
Then the Hamiltonian < f27b is replaced by the slave-boson Hamiltonian as 

Hhc -> Hsb = ~ ^e r e+S r , +^V r b+b r . (166) 

<r,r') r 

A hopping process can be understood as a swapping of an occupied site and an empty site. The occupation 
number operator of site r is b+b r . It should be noticed that the external potential acts only on the particles 
but not on the empty sites. To assure that a lattice site r is either empty or occupied by a boson, we impose 
the constraint 

b+h + e+e r = 1. (167) 

A similar theory for the Bose-Hubbard model has been established in refs. ITUl ITD . In this case, an 
infinite number of operators (b™ ) + , 6" for each occupation number a has to be introduced at each lattice 
site, because multiple occupation is possible. In this respect, the slave-boson approach for hard-core bosons 
is much simpler. However, the hard-core boson model describes a projection of the full Bose-Hubbard 
model to n and n + 1 bosons per site, as discussed in the Introduction ! 1.41 

The grand canonical partition function of the system can be expressed as a functional integral with two 
complex fields b r {r) and e r (r). For the following mean-field calculation, we use the classical approxima- 
tion here, which only takes into account thermal fluctuations but not quantum fluctuations. This means that 
for the fields in Matsubara representation 



br(r) = -j= 6 r ,„ e lw " r ; e r (r) = -= ^ e r , n e 

with bosonic Matsubara frequencies co n , only the terms with uj = are taken into account, if one assumes 
that 

e r ,w„ ~ e r ^ n ~ , if n ^ 0. (168) 

In other words, the time dependence of the fields is neglected. This is justified if we can assume that 
quantum fluctuations (which are neglected in the classical approximation) are small. 
The constraint |6 r | 2 + |e r | 2 = 1 is enforced by a J-function in the integration measure: 



Z sh = J e- A ^ b *- e ^V[b,b*,e,e*} , 



(169) 
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with 

V[b, b*,e, e*] = (\br\ 2 + |e r | 2 - 1) d& r d&Jde r de* 

r 

and the action 



(170) 



A[b, b*,e, e*} =p\-J2 - Y d E b 

[ r ( r>r ') 



r e r e*,& r < 



(171) 



Here, we consider a space-dependent chemical potential /x r = jj, — V r . 



5.2.2 Two-fluid theory in classical approximation 

The hopping term of the action is of fourth order in the field variables. Therefore it is not possible to 
perform the integration directly. However, it is possible to decouple the hopping term by introducing two 
new fields, a complex field $ and a real field ip, and perform a Hubbard-Stratonovich transformation. The 
fields b and e can be integrated out then, and a mean-field approximation can be applied to the fields <f> and 

The idea of the Hubbard-Stratonovich decoupling is similar to the one used in the previous chapter to 
decouple the fourth order terms of the Grassmann fields. We insert the identity 



const, x e 



-A[b,b* ,e,e*} 



2sip r - 



,(e r ,6 r ) 
with the integration measure 

x>[**,*,d = n 



exp 



s<I> r 

-Mr 



■0 



E*: 



d<&*d<i> r d93 r 

(2tt) 3 / 2 



(172) 



(173) 



Here, J is the hopping matrix ( 176*1 ). The constant factor is of no physical relevance. Like for the paired- 
fermion model which was discussed before, the parameter s takes care of the convergence of the Gaussian 
integral. It has the unit of an energy and should not be too small compared to J. Although the exact identity 
does not depend on s, we will see subsequently that the mean-field equation we will derive, does. This is 
a difference to the previously discussed model, where the result which was derived on the mean-field level 
and on the level of Gaussian fluctuations, did not depend on the free parameter s. 

After substituting the identity ( 1172b into the functional integral ( 1171b . the fields b and e are only of 
second order and can be integrated out exactly together with the constraint. This is shown in Appendix IF. 1 1 
The result for the partition function is 



(174) 



with the new action 



s-J 



(175) 
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and the function 



sinh 



dip r 



s 2 l$ r 



(176) 



Note that the action v4(<I>*, $) does not depend on the real field <p explicitly, because it appears inside the 
function Z' only as an integration variable. 

The form (1174b of the grand canonical partition function can be understood as a two-fluid theory. It is 
shown in Appendices IF.2I and IF. 3 1 that the condensate density is related to the field $ and is given by the 
relation 



n 



lim (* r $;,) , 



r— r — >oa 



(177) 



and that the total particle density at site r is related to the field ip by means of the expectation value 

n r = fo.) + i . (178) 



5.2.3 Mean-field theory 

A mean-field solution is found by minimising the action via the variational principle 5 A = 0, which leads 
to a saddle-point approximation, as it was done for the paired-fermion model. Since the field ip can be 
integrated out (e.g. numerically) inside the function Z' r given in Eq. (1176b , minimization has to be done 
with respect to the complex field $ only: 



dA 



dA 



= . 



This yields the mean-field equation 



E 



j 



8 



d(|$ r p) 



$ r = 



(179) 



(180) 



In the case of a spatially constant field without external trapping potential, i.e. if we assume that $ r = <i>o 
and /i r = fj,, the mean-field equation is 



1 



8 



\ogZ' = . 



s + J /?d(|$o| 2 ) 

If the field $ is varying only very slowly between neighbouring lattice sites, we can approximate 



E 



j 



J 



(s + J) 



2 Z-^i 



J5 r 



J v 



(181) 



(182) 



In Fig. Q~2]the phase boundary between the BEC and the non-condensed phase is plotted for different values 
of s. The phase boundary solves Eq. (11811 ) for $o = 0, and has been calculated numerically. 

One can see that the BEC phase forms a "bubble" in the phase diagram, if s/J > 1. This behaviour 
is unexpected because the BEC phase should become narrower, if temperature is increased. This means 
that for too large values of s/J the mean-field theory seems to be incorrect. However, it turns out that the 
absolute minimum of the action with respect to s at constant J, \i and (3 occurs at values of s/J < 1. 
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Fig. 12 Phase boundary between the BEC and the non-condensed phase for s/J — 3 (long dashes), s/J = 1 (short 
dashes), s/J = 0.2 (solid line). Compare these graphs with the graph on the right hand side of Fig. [9] where the 
critical temperature of the mean-field result for the paired-fermion model is plotted. 
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Fig. 13 Total particle density and condensate density for zero temperature (thick lines) and for non-zero temperature 
(thin lines, s/J = 1/5.5) against chemical potential |63|. Compare this graph with the result for the paired-fermion 
model plotted in Fig. 1101 



It is possible to find an exact solution for zero temperature, which does not depend on s. This calculation 
is shown in Appendix |F.4| Two phase boundaries are found: A boundary between the BEC and an empty 
phase with fj, c = — J and a phase boundary between the BEC and the Mott insulator with /i c = J. It is 
identical to the zero temperature mean-field result in Eqs. d 1541 ) and J155I ) that was found for the paired- 
fermion model, and agrees with it qualitatively at finite temperatures (see Fig. [13}. When temperature 
increases, results strongly depend on s. 
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5.2.4 Quasiparticle spectrum 

We get the quasiparticle spectrum from the Gaussian fluctuations, the same way as it was done for the 
paired-fermion model. We write 

$ r = $ + <5$ r , $; = % + 5<s>; , 

and assume that the fluctuations SQ, <5$* about the mean-field solution $o are small. Substituting this 
expression into the action ( 1175b . and expanding it up to second order in the fluctuations, one finds 

A = P^fj\*o\ 2 -logZ'(\$ \ 2 )- f X> $ "«S)£2 ( HI ) , d83) 



r,r 



with the matrix 



Here, we have introduced the abbreviations 

h2 := -ldW) logZ ' 



S + J 



(185) 
(186) 



$=3>o 



and used the approximation in Eq. (11821 ). The matrix Q has no time-structure because of the classical 
approximation. To find the Green's function of quasiparticles, we artificially introduce the imaginary time 
by writing 



rr ' V $ 2 ~a 4 5 rr , J6 rr ,+j„,-h& +a2 + | $0 |2 Mrr/ y ' 

in analogy with the Bogoliubov theory. After a Fourier transformation it leads to the Green's function 



G {k >" n) -VTW{ i^„ e k+ (-a!(a 2+ 2 W ) J ' (188) 



which is equivalent to the matrix (II 14b . and et is the lattice dispersion ( I58l l. The quasiparticle spectrum is 
given by the poles of Q, and can be found by performing the analytic continuation ifkj n — ► and solving 
the equation det Q^ 1 =0. We find solutions for both the BEC phase and the non-condensed phase: 

In the BEC phase, where l^ol 2 > 0, the coefficient a,2 vanishes, because $o solves the mean-field 
equation (118 It . which is equivalent to 5,2 = 0. The solution is 



^k = W£k[2 ( ' S+ / )2 «4|1>o| 2 +ek) ■ (ISO) 



It is gapless and agrees with the Bogoliubov spectrum dl 16t , when we identify the condensate density with 
no = s 2 |<I>o| 2 /(s + J) 2 , and the interaction constant with g = (s + J^a*/ s A . The coefficient 0,4 depends 
on both temperature and chemical potential. Its zero-temperature result is given in Eq. J247t of Appendix 
IF.4I In the dilute gas (i.e. near the phase transition to the empty phase) where uq <C 1, we find at zero 
temperature for the interaction constant the result g w 2 J. 
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In the non-condensed phase, where \&o\ 2 = and a-2 7^ 0, the quasiparticle spectrum is gapped, in 
agreement with the findings of the paired-fermion model: 

E k = e k + A , (190) 

with the gap A = (s + J) 2 a2/s 2 . At zero temperature and near the phase transitions, we find the result 
A = \fi— fi c \+0((n — He) 2 ) which is identical to the zero-temperature result ( 11561 ) for the paired-fermion 
model. 



5.2.5 Renormalized Gross-Pitaevskii equation 

In this section we will derive a mean-field equation which is appropriate to describe the BEC as well as the 
Mott insulator in a strongly interacting Bose gas, and which is similar to the stationary Gross-Pitaevskii 
equation. The mean-field equation for a hard-core Bose gas in an optical lattice within the slave-boson 
approach is given by 



(s + jy 



4>r 



$ r 

s + J (3 



d 



d(\t> r \ 2 ) 



iogz; 



$ r = 



(191) 



This we get by applying the approximation J 1 82b in Eq. (11 80b . However, it also possible to describe a 
system of strongly interacting bosons without lattice potential within this approximation. Therefore we 
perform a continuum approximation of the hopping term: If the lattice constant a is so small that the order 
parameter $ r varies only slowly over neighbouring lattice sites, we can treat the 3-dimensional lattice 
approximately as a continuum: 



JS r 



4\. 



Ja 2 A $ r 
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j'=i 



2$ r 



Ja 2 



V 2 $ r 



When working on the continuum, we rescale the order parameter by 

$(r) := a~ 3/2 $ r , 
such that the action d 175b can be written as 



(192) 



(193) 
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The order parameter is normalied to the number of condensed particles by 

„2 



No 



(s + jy 



J |$(r)| 2 d 3 r . 



(194) 



(195) 



The replacement d 193b has also to be made inside the function Z', of course. The corresponding mean-field 
equation for the continuum is 



6 (3s 2 9(a 3 |$(r)| 2 ) 



logZ'(r) 



$(r) = 



(196) 



The parameters can be identified with those of the conventional GP equation: The mass m of the particles 
is given by the hopping constant J and the original lattice constant a via 
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Fig. 14 Coefficients fin and of the renormalized GP theory plotted against the chemical potential fi. All parameters 
are normalised by the inverse temperature 0. The tunneling rate was chosen to be PJ = 5.5 and the free parameter 
was chosen as s = IcbT. 



In the continuum a looses its identity as lattice constant, but describes a characteristic length scale that can 
be interpreted as the spacial extension of a boson. Thus, it should be of the same order of magnitude as the 
s-wave scattering length a s . 

If the order parameter $ is small, we can expand the potential part of the action up to fourth order: 



(s + J)a 3 |$(r) 
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Z'(r)e- 
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where we have introduced the coefficients 
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They depend on /i, J, /3, and |$(r)| 2 . Further, we introduce the rescaled order parameter 



*R(r) 



s + J 



$(r) 



With these coefficients, the full mean-field equation ( 1196b can be approximated by the equation 

Jo 2 



-^j-V 2 - fi R + g R \^ R (r)\ 2 



*n(r) = 



(198) 

(199) 
(200) 
(201) 

(202) 
(203) 



This equation has the same form as the conventional stationary GP equation, where fi R and g R play 
the role of a renormalised chemical potential and a renormalised interaction constant, respectively. Their 
dependence on /i is shown in Fig. [14] Therefore we refer to this equation as a "renormalised GP equation" 
[72|. The zero temperature limits of the coefficients are calculated in Appendix IF.4I see Eq. ( 12481 ). Near 
the phase transition to the empty phase, i.e. in the dilute regime, where \i = — J + A/x, A/j, <C J, we find 
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ur = A/i + 0(Afi 2 ). Thus, in the limiting case of a dilute BEC and zero temperature, the renormalised 
GP equation goes over to the conventional GP equation with the interaction parameter g = g^ = 2a 3 J. 
While <?r is always positive, ur can change sign. A BEC exists if jUr > 0, otherwise the order parameter 
vanishes. The phase transition between the BEC and the non-condensate phase is given by the relation 
«r = 0, which is equivalent to Eq. ( 11811 ) in a translational-invariant system. Inside the BEC phase, 
/xr increases linearly with increasing «, reaches a maximum and decreases again until the condensate is 
destroyed totally due to strong interaction effects. 

6 Discussion 

6. 1 Comparison of the results 

The main results that we found for the one-dimensional model, the paired-fermion model, and the slave- 
boson model, will be summarized and discussed in this section. All three models give more or less the same 
physics at zero temperature, with an empty phase, a phase with a particle number per lattice site between 
and 1, and a Mott insulator. Their common features and differences shall be pointed out in detail. 

6.1.1 Phase diagram, total density and condensate density 

At zero temperature, the exact solution of the one-dimensional model exhibits three phases in the trans- 
lational invariant case, as shown in Fig. |6]in the J-/i plane: An empty phase which contains no particles 
in equilibrium (physically speaking, it costs energy to put a particle into the system), an incommensurate 
phase with a particle number per lattice site ?i to t between and 1, and a Mott-insulator with n to t = 1- The 
same zero-temperature phase diagram has been found for the paired-fermion model (see picture (a) in Fig. 
|9]l and the slave-boson model on the mean-field level. The only difference is that for the three-dimensional 
models, the incommensurate phase is a BEC, whereas in the case of the one-dimensional model there is 
no BEC but only a long range correlated phase. This is a consequence of the Mermin-Wagner theorem 
[8 9|. At non-zero temperatures, the empty phase and the MI are affected by thermal fluctuations, and 
they have no clear phase boundary any more. However, the three-dimensional systems still have a single 
phase boundary between a BEC with a non-zero order parameter, and a non-condensed phase where the 
order parameter vanishes. The shape of this phase boundary depends on temperature (see picture on the 
right hand side of Fig. |9]for the paired-fermion model). 

For the one dimensional model, the total particle density at T = and T > is shown in Fig. [5] At 
T = 0, the derivative dn to t / dfj, diverges at the phase transitions between the BEC and the empty phase 
and the BEC and the MI phase. The sharp transitions are "washed out" at finite temperatures. 

The zero temperature mean-field results for the total particle density and the condensate density of the 
paired-fermion model and the slave-boson model agree with each other and are given in the Eqs. (1154b and 
( 1155b . We find a total particle density which increases linearly with /i. In the dilute regime the condensate 
density is given by uq = n to t — C(^tot)- If we neglect the terms of order n^ ot , this is in agreement with 
Gross-Pitaevskii theory which assumes that all particles are condensed in this regime. In the absence of a 
trapping potential, a solution of the stationary GP equation is given by 

no = - ■ (204) 
9 

This describes a linearly increasing condensate density no with respect to the chemical potential. Although 
it takes the repulsion into account by a factor 1 j g which is decreasing with increasing interaction constant 
g, the saturation of no cannot be seen in this solution. From the physical point of view, in a realistic de- 
scription for large densities, the particle density must saturate because there is a finite scattering volume 
around each particle. Furthermore, for increasing particle density, the condensate density should reach a 
maximum and for even larger densities, decrease again until its total destruction, because of the increasing 
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interparticle interaction. This is the behaviour that we found for the slave-boson and the paired-fermion 
model in mean-field approximation. A similar behaviour has also been found by variational perturbation 
theory ll73ll . and diffusion Monte Carlo calculations ll74l . In order to describe condensates at higher densi- 
ties, the second order term in the low-density expansion of the energy density has been taken into account 
which leads to a modified GP theory 1141 174117511761 . 

At non-zero temperatures the phase boundaries of the empty phase and the MI are not well defined any 
more, like in the one-dimensional case. The region of BEC shrinks and the condensate density decreases. 
Non-zero temperature results of the paired-fermion model and the slave-boson model are very similar but 
not identical (compare the figs. [TQland[T3l>. This is a consequence of the different mean-field approaches. 
The effect of quantum fluctuations on the zero-temperature result has been studied for the paired-fermion 
model. A condensate depletion was found, but the critical points were not affected (see Fig. [TTJ. 

6.1.2 Excitation spectrum 

The spectrum of quasiparticle excitations is found on the level of Gaussian fluctuations. For the paired- 
fermion model, and the slave-boson model, the expressions for the quasiparticle spectra are summarised 
in the subsequent table: 





in the BEC phase 


in the non-condensed phases 








paired-fermion model 






£k + mI - J 


slave-boson model 




^2^a 4 |a> | 2 + 6 k ) 


£k + (s + J) 2 a 2 /s 2 



Here, ek is the free-particle dispersion relation in the optical lattice, given by Eq. (pgD , We find a spectrum 
which is linear for small wave vectors k in the BEC phase, whereas the spectrum has a gap in the non- 
condensed phases. The gapless spectrum in the BEC phase is caused by a Goldstone mode due to a 
broken global U(l) symmetry lfl5l . The result given for the paired-fermion model is only valid at zero 
temperature. The gapped spectrum is found both in the empty phase and in the MI phase. The result for 
the slave-boson model depends implicitly on temperature via the coefficients a 2 and a 4 given in Eqs. d 1 85b 
and ( 11861 *, and it also depends on the non-physical parameter s. 

We have shown that the zero-temperature results of all three models inside the BEC phase and near the 
phase boundary to the empty phase (/i + J <C J), agree with the Bogoliubov result 

Ek = \J e k (2^i + e k ) . 

The only difference is that the chemical potential is shifted (p, — > /j, + J), because the phase transition in 
Bogoliubov theory is given by [i = instead of fi = — J for the two three-dimensional models. The region 
near the phase transition to the empty phase is the weakly interacting regime, therefore Bogoliubov theory 
is applicable there. The interaction constant was identified as g = 2a 3 J (where the lattice constant a was 
set to 1 in the lattice models). 

The gapped spectrum in the MI that was found in the paired-fermion and slave-boson models is of the 
form 

E k = e k + A . (205) 

We have shown that in the MI phase, near the phase transition to the BEC phase, the gap is given by 
A = fi - J. 

For the one-dimensional system, the excitation spectrum in the incommensurate phase can be found 
indirectly by means of the Feynman relation and is given in Eq. (1991 . It is linear for small wave-vectors k, 
like in the BEC phase of the three-dimensional systems discussed above. 
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6.1.3 Static structure factor 

The static structure factor is defined as the Fourier transform of the equal-time density-density CF, as it 
is defined in Eq. (l54l . At zero temperature it is related to the quasiparticle excitation spectrum via the 
Feynman relation 



where the identification f? /2m = Ja 2 /2d can be considered for a lattice system (in this case m = rn* 
is the band mass as defined in Eq. (|59l>). For the weakly interacting Bogoliubov gas the density-density 
CF was calculated explicitly on the level of a Gaussian approximation. It shows an algebraic decay with 
1/r +1 , where d is the dimension. The result for the static structure factor agrees with the Feynman 
relation. For the one-dimensional system the density-density CF, and therefore the static structure factor, 
were calculated exactly in the incommensurate phase, and agree with results from the literature. In the MI 
phase it vanishes. 

6.2 Comparison with results from the Bose-Hubbard model 

In previous calculations, performed on the Bose-Hubbard model, each phase requires its own specific 
mean-field approach ll55l 1771 or a single one close to the phase boundary [35|. Within a Bogoliubov 
approximation to the Bose-Hubbard model the quasiparticle spectrum in the BEC phase was found as 



where U is the interaction parameter and no is the condensate density. In contrast to this expression, we 
found for the spectrum the expressions in the table in section 16.1.21 These expressions do not agree in the 
limit U — > oo. Thus our hard-core Bose gas cannot be described within the Bogoliubov approximation 
to the Bose-Hubbard model by simply sending U to infinity. On the other hand, our results are in good 
agreement with a variational Schwinger-boson mean-field approach to the Bose-Hubbard model, which 
describe the phases near the phase transition, by sending U to infinity (35). In the large-[/ limit of the 
Bose-Hubbard model, multiple occupation of lattice sites is prohibited because it cost a large amount of 
energy. Therefore one can assume that in this case, the bosons behave like hard-core bosons. 

The results for the excitation spectrum in the Mott-insulating phase from the paired-fermion model and 
the slave-boson model are consistent with the spectrum that was found for the Bose-Hubbard model in 
the large-C/ limit. Inside the first Mott lobe, which is the equivalent to the MI with filling n to t = 1 for 
hard-core bosons, the latter is given by the expression [55 70"l[35l 



which describes two branches: One ("+" sign) is assigned to quasiparticles and one ("— " sign) to quasi- 
holes. It depends on the interaction parameter U. For our hard-core bosons, only the quasihole branch 
can exist, because the hard-core condition prohibits multiple occupation of lattice sites, in contrary to the 
Bose-Hubbard model, where multiple occupation is possible and allows the creation of particle-hole pairs. 
For large values of U the square root term can be written as 





(206) 




such that we find for the two branches the large-C/ results 




(207) 



(208) 
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The gap of the quasiparticle branch is of the order of U, and in the U — ► oo limit it goes to infinity, 
because the energy to occupy a site with two particles is infinitely large. On the other hand, the terms 
which are proportional to U cancel for the quasihole branch, and its U — ► oo limit is identical to the result 
given in Eq. ( 1205b . Particle-hole excitations cannot be created for hard-core bosons, so the creation of an 
elementary excitation is associated to removing a particle out of the Mott-insulator. This is possible in the 
grand-canonical ensemble, where only the average number is fixed but the number of particles fluctuates. 
Inside the empty phase, the same quasiparticle spectrum was found as for the Mott-insulator, due to the 
particle-hole symmetry. Here, the creation of an excitation is interpreted by putting an additional particle 
into the system. 



7 Conclusion 

In this review, the many-particle problem of strongly interaction bosons in a lattice potential was inves- 
tigated. This is motivated by recent experiments on Bose-Einstein condensates in optical lattices which 
showed the phase transition from a BEC to a Mott-insulator. Three different models are discussed, which 
allow the calculation of the phase diagram, and experimentally observable physical quantities like the total 
density, the condensate density, the quasiparticle spectrum, and the static structure factor. All these models 
have in common that they simulate a strong repulsive interaction by imposing a hard-core condition on the 
bosons, which prohibits a multiple occupation of lattice sites. They are defined by means of the functional 
integral method. 

The first model is a special construction which describes non-interacting impenetrable fermions in a 
one-dimensional lattice. We exploited the well-known fact that such a fermionic system is equivalent to 
impenetrable bosons in one dimension, and that the static structure factors of the fermionic and the bosonic 
system are identical. As the fermions are non-interacting, the model can be integrated out exactly. We cal- 
culated the local particle density, the density-density correlation function and the static structure factor in a 
translational invariant system as well as in a system with a harmonic trap potential. In the translational in- 
variant case, the static structure factor, which is experimentally accessible in Bragg scattering experiments, 
increases linearly for small wave vectors, until it reaches unity and remains constant. The density-density 
correlation function shows characteristic oscillations and decays like 1/r 2 . 

The other two models were applied on a Bose gas in a three dimensional lattice. They were treated in 
mean-field theory. The first one, which was called the paired-fermion model, was constructed by a field of 
pairs of Grassmann variables. It can be seen as an interacting fermionic model. The second one was based 
on a slave-boson approach. A Hubbard-Stratonovich transformation allows to integrate out the original 
fields in both models. This transformation leads to new fields, which are connected to the condensate order 
parameter. A saddle-point approximation provides both a mean-field solution and Gaussian fluctuations. 
The latter contain the information about quasiparticle excitations. For a three-dimensional lattice, the total 
particle density and the condensate density can be calculated in mean-field theory, and the quasiparticle 
spectrum and the static structure factor was calculated on the level of Gaussian fluctuations. The saddle 
point approximations of the two models lead to qualitatively the same results. 

Our results for the one-dimensional model, the paired-fermion model, and the slave-boson model, show 
a particle hole symmetry. At zero temperature, they have a common phase diagram, with one phase bound- 
ary between the empty phase and the incommensurate phase, and one between the incommensurate phase 
and the Mott-insulating phase. If the temperature is non-zero, there is no clear phase transition between 
the empty phase and the Mott-insulator due to thermal fluctuations. While there is no Bose-Einstein con- 
densation in the one-dimensional system, the incommensurate phase is a BEC in the paired-fermion and 
slave-boson model in three dimensions. For the latter two models, the mean-field results for the total 
density and the condensate density agree exactly at zero temperature, at higher temperature they agree 
qualitatively. It was shown that they lead to the Gross-Pitaevskii result in the limit of low temperature, if 
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the density is small compared to the lattice constant. At higher temperatures, we have shown that the slave- 
boson model leads to a renormalised Gross-Pitaevskii equation with temperature dependent coefficients. A 
similar theory could in principle be derived on the mean-field level from the paired-fermion model as well. 
It could be compared to the renormalised Gross-Pitaevskii theory which was derived from the slave-boson 
model. 

The quasiparticle spectra which were found for both three-dimensional models, are gapless (Goldstone 
mode) in the BEC phase. In the dilute regime, they agree with the well-known Bogoliubov result. In the 
empty phase and the Mott-insulator, the quasiparticle spectrum is gapped. Our results agree with results 
which were derived for the Bose-Hubbard model, if the on-site interaction constant U is very large. The 
Goldstone mode in the BEC phase of the paired-fermion model was found as the quasiparticle pole of only 
one eigenvalue of the 4x4 quasiparticle Green's function. Additional massive modes may be found from 
the remaining eigenvalues. 

At zero temperature, the elementary excitations are connected to the static structure factor via the Feyn- 
man relation. In the empty phase and the Mott-insulator, the static structure factor vanishes because of the 
absence of density fluctuations. 

A Finite sums and products 

A.l Bosonic sum 

For bosonic systems, which have a periodic structure in the imaginary time variable, we have to perform 
sums of the type 

W 1 , ^TT"" 1 



This sum is performed by finding the common denominator, which is given by 1 — a M . The numerator 
then is 
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numerator = ^ e ~T lml J^J ( 



1 - ae + T fe 



n=l k^n 



where 



Therefore we find 

M M M 



numerator = e ~^" m a 1 1 e^ L< - 1 1)n = a 1 x e 



n(m — Z+l) 

n— 1 1=1 n,l—l 

M oo 

l — l j=—oo 

With the restriction m = — [M — 1), . . . , M — 1 the "enhanced" Kronecker symbol 6' contributes for the 
two cases 

I = m + 1 if m > 

I = M + m + 1 if m < . 

Finally, this leads to the components of the inverse matrix: 

M 



^ 1 e -— nm 1 f a m if m>0 

^Afi- e¥» " l-a M X I 



M - ) a M+m tf m<Q ■ (209) 
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A. 2 Fermionic sum 

For fermionic systems, which have an anti-periodic structure in the imaginary time variable, we have to 
perform sums of the type 



1 e m 



M 



1 



(n-i) 1 + a 



M 



,M4 



if m > 
if m < 



(210) 



This sum differs from the sum given in Eq. ( 12091 ) only by the substitution a — > ae mm / M and a multipli- 
cation by the factor e mm / M , so the result can be verified easily. 



A.3 Sums with cosines 
The following two sums require the condition \b\ > 1 

M , 



E 



i 



1 (b - % /6 2_ I) M + (b + VW~1) M + 2 



{Mcos(^n)-b (b-Vb^l) M - (b + Vb^T) 



1 cos(ffn) 



i (b-Vb^T) M l + {b + Vtf=T) lw ' + 26 



it/ 



M-l 



E 



-Mco S (^n)-6 (b-VV^T) M -{b + VV^T) 

To perform these two sums the following identities were used: 



M 



2a 2 



cos 



1 



1 1 



cos(x) 



cos(x) 



2a 



a 2 - 1 



12: — a a a e lz: — 1 
11 11 



ae lx — 1 a e lx — a a e 



a ae 



All separate terms can be traced back to the sum given in Eq. ( 12091 l. 



(211) 



(212) 



A.4 Sum for C(k) in Eq. 
We perform the sum 



A I 



C(k) = lim V 



1 



27T1 1 771 



e si" 



M— >OC ^— ' M I AJLLI — 

i=i i e *f — e m 



(e-'-eS (l-^)) - e V'rf 



,2 fe 



Make the following substitutions: 



1 _ ' 61 



/(*) 



b = — Je2M cos — 
M 2 



z + a 



(z + a) 2 — fo 2 z 



With these definitions, the sum is given as 

C(k) = - lim V — / f 



M 



1=1 



(213) 
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The roots of the denominator of f(z) are 



fo2 ^ 

z ± = o± - \/b 2 - 4a . 

2 2 

We perform an expansion into partial fraction and find 

A B (A+ B)z- (Az~ + Bz+) 



z — z T z — z 



with 



(z — z + ){z — z ) 



1 



; B = -- 



2 2y/b 2 - 4a 2 2y/b 2 - 4a 

To perform the sum, we use the following identity which can be traced back to Eq. ( [209) : 



M 

E 

i=i 



i i 



M e^iri - z ± 



l-(^) 



M 



M 

£ e*/(e»') 



, M 



The limit M — > oo can now be performed, by the help of the identities 



lim (z±) M = e" Km ( 1 + ( ± Jcos ~ - /x ) £ + O ( ^ 



M- 



M- 



M 



1 



lim z — 1 : lim A, Z? 

Af— >oo M^oo 2 



M 2 
1 



A/ 



_ _ a(±jcos 



The result is given in Eq. (1871 1. 



A.5 Sum for G in Eq. Qg3) 
We perform the sum 



A/ 



i(i<^o + Xo) 



G = -E 

M «=! 1 + (i**> + Xo)(V5 + XS) - 2e-T(«-l) + e^M) 



(214) 



We define 



a:=l + (i^ + X )(i^+X*), &:=1 -(^ 



/(*) = 



1 



2,2 + oz 2 ' 



The roots of the denominator of f(z) are 

z ± 

An expansion into partial fraction leads to 

1 



1 



= - ( 1 ± Vl - ab 



m = a 



z — z T z — z 



where A = 



1 

2-s/l - ab 



49 



To perform the sum, we use the following identity which can be traced back to Eq. ( 1210b : 



1=1 



^ i 



, M 



M 



A. 6 Product to calculate the determinant of Eq. ( 1120b 
We want to perform a product of the type 



M 



n * 



2tt 



n , |6| > 1 



This can be verified to be equal to 



n 

n=l 

such that the identity 

M 



i (b + Vb 2 - l) (l - (b - Vb 2 - l) e 1 ^") (l - (& - V^ 2 - l) e _i * ' 



1 - aeT™ = 1 - a 



n 

n=l 

can be applied. As a result we find 



. M 



(215) 



A ( 5_C0S (S")) =2- M ^(&+V6 J_ T) M +(6- V / fo JT T) M -2^ . (216) 

B Coherent states for bosons and fermions 

The functional integral representation for bosonic and fermionic systems is constructed of coherent states 
[36|. We denote bosonic operators by a+, a a , and the fermionic operators by c+, c a . The commutation 
relations are 



[a a ,a a ,\_ = 5 aa > , (217) 
[c Q ,c+] + = 6a„r. (218) 
The vacuum state, i.e. the state containing no particle, we call |0). We define coherent states for 
• bosons by means of complex field variables , tj> a : 



3 E K&Z 



10), (01 



(219) 



• fermions by means of conjugate Grassmann variables ip a , ip a , where we require, that the Grassmann 
variables anticommute with the fermionic operators: 

|V») = e-S.^i| )=JI(l-^)|0>, 

a 

(<tl>\ = (0| ^ = (0| 11(1 + ■ (220) 
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For the construction of the coherent state functional integral, the following properties are relevant. They 
can be checked by using the previous definitions and the integration properties of complex, Grassmannian 
and nilpotent variables: 

• Coherent states are eigenvalues of annihilation operators: 

x a \0 = U0, {t\£Z = (Z\L, (22i) 

where x = a, £ = <fi, £ = <fr* for bosons, and x = c, £ = £ = <fi for fermions. 

• Scalar product, where the operator X is built of bosonic, fermionic, or hard-core operators, respec- 
tively: 

(Z\X(x + ,x)\0 = e^^'«X(Z a ,0 , (222) 
where x, £, £ have to be chosen as mentioned above. 

• Closure relation (the unity operator is denoted by 1): 



\<f>) (m^r (223) 

a 

\4>) (4>\Y[d$ a d^ a . (224) 



• Trace of an operator X: 



TrX(ai,a a ) = f e £. O- { ^) J[ (225) 



TrX(c+c a ) = f e-^^(-^\X\i>)Y[d4, a d^ a . (226) 

J a 

Using these identities, the functional integral of the grand canonical partition function 

Z = Tr e -P(Hdi,£ a )-^N(xi,x a )) 

with the Hamiltonian H is constructed in the following manner: We apply the relation for the trace and 
insert the closure relation M — 1 times. Introducing the discrete-imaginary-time index n = 1, . . . , M we 
have 

Z = Je^M^ia^e^^M) fl \in-x) U , (227) 



a ,n 



where a = +1 for bosons and —1 for fermions, and M = 2iri for bosons and 1 for fermions. The minus 
sign inside the scalar product in the fermionic trace gives rise to the anti-periodicity of the fermionic field 
variables. The different sign in the exponent of the hard-core bosonic trace is the reason that the diagonal 
term in the action for hard-core bosons is different from bosonic and fermionic actions. 

The operator in the exponent H(x~^ , x a ) — fj,N(x+ , x a ) can be replaced by its normal ordered from by 
making an error of the order ((3/M) 2 which vanishes for M — ► oo. Applying the eigenvalue property and 
the product property yields 
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with the action 

M 

M 



n—1 \, a 



M 



•ff(£,n+i.&,»n (229) 



and the boundary condition £ Q i = <72£ a) M+i> £a,i = 0"2£a,M+i- 

C Expectation values and Wick's theorem 

An expectation value of an expression in terms of real/complex/Grassmann variables is defined by means 
of Eq. (HTt . A second order expectation value provides the matrix element of the (inverse) Green's matrix 
Q: 

Real variables: (<l>j4>k) = h^jk 

Complex conjugate variables: {4>*j4>k) = GJ^ (230) 
Conjugate Grassmann variables: ("ipjipk) = Gjk 

Forth order expectation values can be calculated via the application of Wick's theorem lfl4l [36l . It can 
be split into products of second-order expectation values and a sum has to be performed over all possible 
pairings (including a sign for Grassmann variables): 

Real var: (<pj<f> k <f>i<p m ) = {(f>j<pk){<l>l<Pm) + (<Pj<t>l)(<f>k<l>m) + (<t>j<Pm)(<l>k<Pl) 

C.conj. var: (cpjcgMm) = <^m><<^> + <^><<^>m> (231) 

Conj. Gr. var.: (i/jjipk^l^m) = (ipjipm) (tpk^l) ~ (">Pjipl) (^fcV'm) 

D Correlations 

The decay of the density-density CF given in Eq. ( 11291 ) is investigated in d — 1, 2, 3 dimensions. For 
convenience we write c :— y / 2(^7+J). We use a cut-off at |q| = Q for the integrals. 

• One dimension: 

r Q iffi 2 r Qr i 

D{r) = J_ Q f e-' dq = — J o q> cos(q>)d q > ~ 

The anti-symmetrical part which is ~ sin(g') does not contribute. 

• Two dimensions with polar coordinates (q, 4>): 



D{r)=l dqq tylj*—*^ ^-3-/ g ' 2 CO S (g')dg'~-3 



Three dimensions with spherical coordinates (q, 6, cj)): 



f Q f 2-K r -l 

D(r)= dqq 2 d</> / d(cos 6) - e lqr cos9 
Jo Jo Jl c 
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E Calculations to the paired-fermion model 

In this Appendix we write out the expression for the Green's function in both cases \<f»\ = and \(f>\ ^ 0. 
Case: \<j>\ = 

Deviation of the effective action due to fluctuations is 



g- 1 



SA e g = 5<f> kiU 5 X k,w ) ( k . n/ , ' i 



(232) 



where 



D{w) = 



1 



1 



\H\—Ui3 " J(3-£k' 

The determinant of the Green's function reads 
dot 0- 1 = - D(u) (^--vr 1 



2 J 



\2J 



(233) 



Case: \(f>\ =f 

Deviation of the effective action due to fluctuations is 



/ 



SA eS = ^2( 6</>w iU ,6xk.,u>,S(l>*_ k> _ u ,8x*-- kt - ul )C 1 



J k,u 



V Sx-k.-u, ) 



(234) 



with the Green's function 



g- l = 



(J) \D{oj) —a ia \ 

iD(ui) jj + D(u) ia a 

—a ia t)^ 1 — D(uj) \D{uj) 

ia a LD(u;) jj + D(u) J 



(235) 



where 



1 ^ 2 + J 2 + 2i M ^ 
^ )= 2' J(J 2 +c 2 ) ' 



1 + J 2 — 2ifxuj 
D{ ~ U) = 2 ' J(J 2 +^ 2 ) ' 



1 



|$| 2 /9 



2 J(J 2 +co 2 )' 
The determinant of the Green's function is 

1 



det^r 1 = 



[2J2(3-e k )] 2 (^ 2 +^ 2 ) 



(236) 
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F. 1 Integration of the constraint 

We perform the integration of the complex fields b and e. The integral factorises such that it can be 
performed for each lattice site r independently. Therefore we will drop the index r here temporarily and 
evaluate the expression 



exp <^-(3sip 2 - (3(e, b) 
The eigenvalues of the 2x2 matrix are 

A± = (3s [<p + \ 



2sip + s s$ 
s$* -fx 



l)de*ded&*d& . (237) 



s 2 m 



A unitary transformation can be applied to the vector (e, b) such that the matrix has diagonal form. This 



does not affect the constraint, because the expression \b\" 



1 remains unchanged after a unitary 



transformation. Therefore the integral is equal to 

de*ded6*d6exp [-(3sip 2 - Ai|e| 2 - A 2 |6| 2 ] 5{\b\ 2 + |e| 2 - 1) 



(2 



it) 2 - / dp p exp [-(3sip 2 - Aip 2 - A 2 (1 - p 2 
2 Jo 



4.7T exp 



1 



2 e -0s V 2 



-\2 



Ai — A2 



0s<? - (3s [ if + £ I + & 



sinh 



/V[fa+5) a + £] 2 + s2 i*i : 



After performing the shift (p + 1/2 — > ip and using the index r again, the integral d237b gives the result 
sinh 



d<p r 



■ s m T 



(238) 



F.2 Condensate density 

In a Bose system in an optical lattice, which is described by a complex field <f) T (t), the condensate density is 
defined by the expression 021 via the concept of off-diagonal long range order. In classical approximation, 
the field does not depend on imaginary time t, and in the slave-boson approach, we replace 



b*e r ; 



e*b r , 



thus we use the definition 

n = lim (6*e x e^,6 X ') . 

x— x '— >oo 

for the condensate density. Here, the expectation value is given by 

(•••) = ±- J ■■ ■ exp[. . .] £>[$*, 3, <p] V[b, b*,e, e*] . 



(239) 



(240) 
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We are interested in the connection between the correlation function (^x'&x') an d t ^ e condensate density. 
For this purpose we integrate out the field $ to transform the correlation function of the field $ back to a 
correlation function of the fields 6 and e. Therefore, we write 



S Or r 1 J p ■ 



for simplicity and perform the integration 



(3 2 s 2 / $ x $ x , exp 



d 



d 



5(& x e x ) d(b x ,e* x ,) 



exp 



(3 +PsJ2 $ r^e r + /3s J2 Ke* r b r 



d 



d 



det — exp 



(3s 2 b*e r v rr /e*,b r 

r,r' 



(3s 2 dot [- 



(3s 2 ^ b*e r e*,b r >v rx v x > r > 



exp 



(3s 2 6*e r w rr 'e*/6 r 



Since we are interested in the limit x — x' — > oo, and the matrix j xx / includes nearest-neighbour hopping 
only, the term w xx / vanishes. This yields for far distant lattice sites x, x' the expression 

($ x $ x >) = s 2 Y^ (b*e r e*,b r ,) v rx v*>r> ■ 



Further we can assume that (b*e r e*,b r >) — (& x e x e*,& X ') for r, x and r', x' nearest neighbours. Using 

s + J 



we get 



and therefore 



lim ($*$ X ' 

x— x'— >oo 



(s + J) 2 , s 



lim (& x e x e*/6 x 



n 



(S + J) 2 X-X'^OG 



x— x — *oo 



lim 



F.3 Total particle density 
The total particle density at site r is given as 

n r = 1 - (M 2 ) , (241) 

where e is the field associated to empty sites. It is possible to express the expectation value of the complex 
field e in terms of an expectation value of the real field ip. To achieve that, let us regard the integration over 
the fields b, e, and <p. After performing the substitution tp + 1/2 — ► cp and dropping the index r, we have 



J dcpe-Hv-iT J P[6,6*,e,e*]|e| 2 exp|-/3(e,6) 



2sp s$ \ / e* 
s$* -/i J I b* 
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Partial integration leads to 
1 

2s~0 

Therefore we find 



dip 



-20s ip 



-ps^-iy J e e *j exp | 



(|e| 2 > = (-U- 



1 



Together with Eq. ( 12411 ) we find for the local total particle density the expression 



(242) 



F.4 Zero temperature limit 

We want to integrate out the function Z' (we drop the index r) given in Eq. d 1 76b for zero temperature, i.e. 
in the limit — > oo. For simplicity we write :— s and perform the limit — > oo instead. Further we 
write a := fi/2s, and x := |(E>| 2 . The function Z' we write as 



Z' = ±*(Z- 

2/3 



Z 



+ ) > 



where 



and 



-0f±( V ,x) 



y/(ip + a) 2 + x 



dip 



f±(f,x) = f 2 ± yj{(p + a) 2 + x 



In the limit (3 — * oo we can calculate the ^-integral Z± exactly by means of a saddle-point integration. 
This is done by expanding the functions f± in second order about their minimum with respect to (p. We 
need partial derivatives 



df±(<P,x) 
dip 

d 2 f±{<fi,x) 
dip 2 



2ip± 
2± 



ip + a 



\/(ip + a) 2 + X 
x 



[(<p + a) 



We determine the extrema of f± : 

df±(ipo,x) 
dip 

which is equivalent to 







y/ (ip + a) 2 + x = =F 



Po + a 
2ip 



x = (ip + a) 2 



1 



1 



(243) 



(244) 



Thus the saddle point approximation for large values of is 

Z± f 



f3\f ± ( Vo , x )+i^l±( Vo , x )( v - Vo y 



y/ (<pq + a) 2 + X 



■ dip 
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-Pf±{<f>a,x) 
(tp a + a) 2 +x [0 d 2 f ± ( Vo ,^ 

" 2 



From Eq. ( 1243b we get 



/±(<A)) 



= 2 



8xf> )¥>o 



2<^ ' ^v? 2 (ipo + a) 3 

where x itself depends on ipo independently via Eq. ( I244l i. For given x there are two solutions for cpo, but 
only the one which is the absolute minimum contributes to Z' for large values of (3. Therefore: 



\ogZ' = log(<^o) - log((^o + a) - g lo S 



d 2 f±(<Po) 
dip 2 



Pf±{fo) + const 



The term that is proportional to (3 dominates all the others, and in the limit /?^oowe find the exact result 

1 



lim -= log Z' = - 

p^oo (3 

1 d 

=> lim — — log Z' = - 

P^oo (3 ox 



f±( t Po) 
d/±(yo) 



d^ 



dx 



The derivative of (po with respect to x we get from Eq. ( 1244b by means of the implicit function theorem: 



d(f _ _ 

da; (ip + a) (4^0 + a ) 



Therefore: 



1 d 
lim loeZ' 



p() 



/5^oo f3 dx <po + a 

Together with the mean-field equation (11811 >. we find the zero temperature result in the condensed phase 
(i.e. where x > 0): 

=> <Ai 



JL 

2J 



s + J (fo + a 

For the order parameter we find from Eq. ( 1244b in the condensed phase: 

2 



1*1 



Thus the condensate density by the definition in Eq. ( 1177b is: 



n 



J? 



|$|2 = 



if-J< M <J 



else 



and because of (ip) — (po the total particle density by the definition ( 1178b is: 



"-tot = <A) + 



if fi < - J 

if -J < [i < J 

if J < n . 



(245) 



(246) 



To determine the coefficient £4 in Eq. ( 1186b . we need the second derivative of log Z with respect to x: 

dipo 
dx ' 



1 d 2 

/§^oo (3 dx 2 ° S 



d<p p-roo (3 ox 



1 



"Wo 



s (ip + a) 3 (4ip 3 + a) 



57 



With the above solution this yields 



1 d 2 s 4 
lim - ■— log Z' = 2 J- -r 

/3^oo f3 dx 2 B (S + J) 4 



1-4 



s + J 



(247) 



With these results we also find the zero temperature expressions for the renormalised coefficients ( 1200b and 



(s + J) 2 

MR = -(« + •/) + ,|i ; 5R = 2a 3 J 



(248) 
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